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PREFACE TO THE FIRST EDITION. 


Of all the celestial bodies whose motions have formed 
the subject of the investigations of astronomers, the Moon 
has always been regarded as that which presents the 
greatest diffici^lties, on account of the number of inequali- 
ties to which it is subject; but the frequent and important 
applications of the results render the Lunar Problem one 
of the highest interest, and we find that it has occupied 
the attention of the most celebrated astronomers from the 
earliest times. 

Newton’s discovery of Universal Gravitation, sugg<»8ted. 
it is supposed, by a rough consideration of the motions of 
the Moon, led him naturally to examine its application to 
a more severe explanation of her disturbances; and his 
Eleventh Section is the first atteinpt at a* theoretical 
investigation of the Lunar inequalities. The results he 
obtained were found to agree very nearly with those 
determined by observation, and afforded a remarkable 
confirmation of the truth of his great principle ; but 
the geometrical methods which he had adopted seem 
inadequate to so complicated a theory, and recourse has 
been had to analysis for a complete ‘ determination of the 
disturbances, and for a knowledge of the true orbit. 
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The following pages will, it is hoped, form a proper 
introduction to more recondite works on the subject: the 
difficulties which a person entering upon this study is most 
likely to stumble at have been dwelt upqn'at considerable 
length, and though different methods of investigation have 
been employed by different astronomers, the difficulties met 
with are nearly the same, and the principle of successive 
approximation is common to all. In the present work, 
the approximation is carried to the second order of small 
quantities, and this, though far from giving accurate values, 
is amply sufficient for the elucidation of the method. 

The differences in the analytical solutions arise from 
the various ways in which the position of the moon may 
be indicated by altering the system of coordinates to which 
it is referred, or again, in the same system, by choosing 
different quantities for independent variables. 

D’Alembert and Clairaut chose for coordinates the pro- 
jection of the radius vector on the plane of the ecliptic and 
the longitude of this projection. To form the differential 
equations, the true longitude was taken for independent 
variable. 

To determine the latitude, they, by analogy to Newton’s 
method, employed the differential variations of the motion 
of the node and of the Inclination of the orbit. 

Laplace, Damoiseau, Plana, and also Hcrschel and Airy, 
in their more elementary works, have found it more con- 
venient to express the variations of the latitude directly, 
by an equation of the same form as that of the radius 
vector. 

Lubbock and Ponlecoulant, taking the same coordinates 
of the Moon’s position, make the time the independent 
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CHAPTER 1. 

INTRODUCTION. 

Before proceeding to the consideration of the moon’s 
motion, it will be desirable to say a few words on the law 
of attractions, and on the peculiar circumstances which enable 
us to simplify the present investigation. 

1. The law of universal gravitation, as laid down by 
Newton, is tliat ‘‘ Every particle in tlte universe attracts every 
other particle^ with a force varying directly as the mass of the 
attracting particle and inversely as the square of the distance 
between themf 

The truth of this law cannot be established by abstract 
reasoning ; but as it is found that the motions of the heavenly 
bodies, calculated on the assumption of its truth, agree more 
and more closely with the observed motions the more strictly 
our calculation^ are performed, we have every reason to 
consider the law as an established truth, and to attribute 
any slight discrepancy between the results of calculation and 
-observation to instrumental errors, to an incomplete analysis, 
or to our ignorance of the existence of same of the dis- 
turbing causes. 

Of the last source of error there is a remarkable ex- 
amnle connected with the discovery of the nlanet Nentune, 
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which became known to us, as one of the bodies of our 
system,* solely by means of the perturbations it produced in 
the calculated places of the planet Uranus. These pertur- 
bations were too great to be attributed wholly to errors of 
instruments or of calculation ; and therefore, either the law 
of universal gravitation ceased to hold for a body so remote 
as Uranus, or else some unknown cause was disturbing the 
path of the planet. The first supposition was too repugnant 
to an astronomer of the nineteenth century to be entertained 
until all others had failed, — and the second supposition led 
to the detection of Neptune. The distinguished names of 
Adams and Le Verrler will be for ever connected with the 
history of this planet, and their solution of the difficult in- 
verse problem: — “Given the perturbations caused by an 
unknown planet, determine, on the assumption of the truth 
of Newton’s law, the orbit and position of this disturbing 
body,” — will always be considered one of the triumphs of 
Mathematical Science. 

Another remarkable instance of the vindication of Newton’s 
law, where it seemed at first to be at fault, will be found 
mentioned in Chap. Vi. in connection with the motion of 
the Lunar Apogee. The error was owing to an incomplete 
analysis; but, when the calculations were carried out more 
fully, the strongest confirmation of the law was afforded by 
the exact agreement of the results with observation. 

We shall proceed to apply the law of gravitation to in- 
vestigate the circumstances of the moon’s motion ; and then 
see how it will enable us to assign her position at any time 
when suflScient data have been obtained by observation. 

2. The problem would be one of extreme, if not insur- 
mountable difficulty, if we had to take into account simul- 

* It had been seen by J)r. Laimmt at Mnnicli, one year before its being known 
to be a planet. ‘‘ Solar System, by J . R. Hind.” 
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taneonsly the actions of the earth, sun, planets, &c., on the 
moon ; but fortunately the planets are so small or so distant 
that their action may be neglected — at any rate to the order 
of our present approximation — and the attraction of the 
earth, on account of its proximity, is very much greater 
than the disturbing^ action of the sun notwithstanding his 
enormous mass. V^c may therefore treat the question as 
that of one body, the moon, revolving about another body, 
the earth, and continually disturbed by a third body, the 
sun. This is the celebrated problem of Three Bodies. 

Stated in this general form its exact solution has hitherto 
defied the powers of the analyst. The disturbing action 
of the sun can be expressed without difficulty in terms of 
the masses and distances, but the integration of the diffe- 
rential equations cannot be effected. If however, instead 
of retaining the whole disturbing force of the sun, which, 
as we have said, is small compared with the earth’s action, 
we expand its expression in a series and neglect all very 
small terms, it is then found possible to obtain a solution. 

This breaking up of the expression for the sun’s action is 
obviously equivalent to a breaking up of the sun’s force; 
and if we afterwards wish to take into account the very 
small forces which have been thus neglected, or the small 
disturbances due to the action of the planets, the following 
principle will shew that we may do so by considering their 
effects separately, — and that the algebraical sum of all the 
disturbances so obtained will, to an order of approximation 
far beyond that which we contemplate, be the same as the 
disturbance due to their simultaneous action. We shall find 


♦ Since the sun attracts both the earth and moon, it is clear that its effects on 
the moon’s motion relatively to the earth, or the disturbing force, will not be tlie 
same as the absolute force on either body. Tlie absolute force of the sun on 
the moon is more than double that of the earth (Art. 99), but the relative force 
does not exceed of the earth’s action (Art. 25). 



4 


LUNAR THEORY. 


an application of this principle in the investigation of the 
parallactic inequality. 

Principle of Superposition of Small Disturbances, 

3. Let a particle be moving under, the action of any 
number of independent forces, 
some of which arc very small, 
and let A be the position of the 
particle at any instant. Let two of these small forces 
be omitted, and suppose the path of the particle under the 
action of the remaining forces to be AP in any given time. 

Let AP^ be the path which would have been described 
in the same time if also had acted; AP^ differing very 
slightly from A P, the disturbance being PP^, 

Similarly, if instead of w, had acted, suppose AP^ to 
represent the disturbed path, being the disturbance, 

{AP^ ylP,, AP,j, are not necessarily in the same plane, nor 
even plane curves). 

Lastly, let A Q be the actual path of the body when both 

and act. Join P^Q^ 

The two disturbances PP^ and P,§, being due to the 
action of the same small force on the slightly different 
paths AP and AP^^ must be very nearly equal both in 
magnitude and in direction. What difference exists between 
them must be a very small fraction of either disturbance, 
and may be neglected compared with the original path. 
We may, in fact, expect that this difference which is the 
disturbance due to of a disturbance due to will be 
of an order compounded of the orders of the two. Thus, 
if the disturbance due to w, and be respectively of the 
second and fourth orders, the difference between P^Q and 
PP^ would probably be of the sixth order compared with 
AP, Therefore Pf} may be considered parallel and equal 
to PP,. 
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Hence the projection of the whole disturbance PQ on any 
straight line, being equal to the algebraical sum of the pro- 
jections of PP^ and P^Q^ will be equal to the algebraical sum 
of the projections of the separate disturbances PZ^. 

jNext, let there be three small disturbing forces 7/1^^ 

We may consider the joint action of the two as one 

small disturbing force ; therefore, by what precedes, the total 
disturbance along any axis will be the sum of the separate 
disturbances of 7 /q and of the system ; but this last is 

the sum of the separate disturbances of 7 n^^ and 7 n^ : therefore 
the whole disturbance equals the sum of the three separate 
disturbances. 


4 . This reasoning can evidently be extended to any 
number of forces; and if a*, y, z be the coordinates of the 
disturbed particle, <f> (a;, y, z) any function of a:, y, 2? ; the 
disturbance produced in (f) (a;, y, z) will be 


20 (a;, ?/, z) = Sx -f Sz^ omitting ( 2 u?)*, &c., 

where 2 a;= 2 a:j + 2a;2+...= sum of disturbances along axis of x 

due to separate forces, 

...= along axis of y, 

221 = 22;^ -f 22:2+... = along axis of 0; 


therefore 8 ^ [x, y, z) = Sx, + Sy, + Sz, 


dy 


= 8 ,^ [x, y, z) + 8 ,^^ [x, y, z) + &c., 

or total disturbanee equals sum of separate disturbances, 
which establishes the principle. 

yince 0 (a;, z) may be the radius vector, the latitude 
or the longitude of the disturbed body, it follows that the 
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total disturbance in any of these elements is the sum of the 
partial disturbances. 

5. But we must still prove another proposition, without 
which the problem, though so far sim^ificd, would scarcely 
be less complicated than in Its most general form. 

Newton’s law refers to 'particles^ whereas the sun, earth, 
and moon are large nearly spherical bodies, and it becomes 
necessary to examine the mutual action of such bodies. 
Now, it happens that, with this law of force, the attraction 
of one sphere on another can be correctly obtained, and 
leaves the question in exactly the same state as if they were 
particles. [Frincip, lib. I. prop. 75.) 

Attractions of Spherical Bodies. 

6. Let P be a particle situated at a distance OP=^a from 
the centre of a sphere of uniform 
density p and radius c. The 
particle being without the sphere, 
a> c. 

Let the whole sphere be divi 
ded into circular laminae by planes 
perpendicular to OP. Let /SO he 
one of these. PS = £c, PQ = and thickness of lamina = Sx* 

Next, let this lamina be divided into concentric rings. 
Let JRS—r be the radius of one of these rings and Sr its 
breadth, z BPO = ff ; 
therefore r = x tan 0, 

Sr = xsec^0,S0. 

The attraction of an element B of this ring on the particle 
P will be along PP, and the resolved part 

of Ihi. along PO wll be co8«. 

X sec 
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But the resultant attraction of the whole ring will clearly 
be the sum of the resolved parts along PO of the attractions 
of its constituent elements ; therefore, 

attraction of ring ~ - cos 6 — 27rp sin Ohx . 80 : 

° X sec 0 ^ ^ ^ 

rcos-i- 

therefore, attraction of whole lamina SQ=27rpSx siu0.J0 

='2TrpBx(l - . 

Again, 2 ;’ = as* + c* — (a — xY = 2aa; — (a* — c*) ; 
therefore zSz — aSx^ 

and attraction of lamina = 27rp 5 

/. attraction of whole sphere = 27rp I ^ ^ „ — — - , ^ 

^ ^ (2a Ga* 2a* j 

(from j 2 ? = a — c to z — a^-c) 

dTTpe^ 

3 a* 

M 




where AT = mass of sphere = 


Air pc 


Hence, the attraction of the whole sphere is precisely the 
same as if the whole mass were condensed into its centre. 

Cor. 1. The attraction of a shell (radius c and thickness 
Sc) will be obtained from the preceding expression by diffe- 
rentiating it with respect to c, and is 

« , 47rpc*.8c mass of shell 

attraction of shell = — ^ — = , 

a* a* ’ 

the same as if the mass were collected at its centre. 
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Cor. 2, Therefore, the attraction of a heterogeneous sphere 
on an external particle loill he the same as if the whole mass 
were condensed into its centre^ provided the density he the same 
at all points equally distant from the centre^ for then the 
whole sphere may bo considered as tne aggregate of an 
infinite number of uniform shells. 

7. Let us now consider the case of one sphere attract- 
ing another. Suppose P in the preceding article to be an 
elementary particle of a sphere M\ whose centre 0* suppose 
at a distance a from 0, Then, since action and reaction 
are equal and opposite, P will attract the whole sphere M 
just as it would do a particle of mass M placed at 0. The 
same is true of all the elementary particles which compose 
the sphere M\ therefore the sphere ilP will attract the sphere 
M as if the whole mass of the latter were condensed into its 
centre 0 ; but the attraction of the sphere JW on a particle 
0 is the same as if the attracting sphere were condensed into 
its centre O' ; therefore, 

Two spheres attract one another as if the whole matter of 
each sphere were collected at its centre. 

8. This remarkable result, which, as may be shewn, holds 
only when the law of attraction is that of the inverse square 
of the distance, or that of the direct distance, or a com- 
bination of these by addition or subtraction, reduces the 
problem of the sun, earth, and moon to that of three par- 
ticles. The slight error due to the bodies not being perfect 
spheres is here neglected, being of an order higher than 
that to which we intend to carry the present investiga- 
tion : this error, however, though very small, is appreciable, 
and when a nearer approximation is required, it becomes 
necessary to have icgard to this circumstance, (See Ap- 
pendix, Art, 109). 
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CHAPTER 11. 

MOTION KELATIVE TO THE EARTH. 

9. When a number of particles are in motion under their 
mutual attractions or other forces, and the motion relatively 
to one of them is required, we must bring that one to rest 
and then keep it at rest, without altering the relative motions 
of the othex’s with respect to it. 

Now, first, the chosen particle will be brought to rest 
by giving it at any instant a velocity equal and opposite to 
that which it has at that instant; secondly, it will bo kept 
at rest by applying to it accelerating forces equal and 
opposite to those which act upon it. 

Therefore give the same velocity and apply the same 
accelerating forces to all the bodies of the system, and t.^ cir 
absolute motions about the chosen body, which is now at rest, 
will be the same as their relative motions previously. 

Problem of Two Bodies. 

10. As the sun disturbs the moon’s motion with respect 
to the earth, it is important to know what that motion would 
have been if no disturbance had existed, or generally : — 

Two bodies attracting one another with forces varying 
directly as the mass and inversely as the square of the dis- 
tance^ to determine the orbit of one relatively to the other. 

Let M' be the masses of the bodies, M' being the body 
whose motion relatively to M is required, r the distance 

c 
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between them at any time and 6 the angle between r and 
some fixed prime radius. 

. M 

Tbe accelerating force of M on M' equals -j- acting towards 

il/' . ’ 

while that of M' on M equals - in the opposite direction. 

Therefore, by the principle above stated, we must apply to 
both M and M' accelerating forces equal and opposite to this 
latter force, and M' will move about M fixed, the accelerating 

force on M being if ^ = M + M ' and r = i . 

r u 


Hence, 


dSi , fM 


d6 

where h = r^ -r— twice the area described In a unit of time, 
dt ^ 


and integrating 

M = +ecos(^ 


e and a being constants to be determined by the circumstances 
of the motion at any given time. 

This is the equation to a conic section referred to its 

/i* 

focus, the eccentricity being e, the semi-latus rectum — , and 

the angle made by the apse line with the prime radius a. 

By observations made on the actual path of the moon, 
which cannot difier very widely from its undisturbed orbit, 
we infer that this latter would be an ellipse with an 
eccentricity about 

In the same manner. If the sun and earth be the two 
bodies considered, the relative orbit would be an ellipse with 
an eccentricity about -ic- 


11. The angle 0 — a behveen the radius vector and the 
apse line is called the tru'^ anomaly. 
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If 71 is the angular velocity of a radius vector whicli, 
moving uniformly, would accomplish its revolution in the 
same time as the true one, both passing through the apse 
at the same instant ; then a is called the ??iea?i 

anomaly where s is a constant depending on the instant 
from which the time is reckoned, its value being the angle 
between the prime iidius and the uniformly revolving one 
when ^ = 0. 

Thus, let MT be the fixed line or prime radius, 

A the apse, 
iir the moving body 

at time 

MfjL the uniformly re- 
volving radius at same time, the 
direction of motion being repre- 
sented by the arrow. 

And, let Ml) be the position of Mfi when t = 0, 

then XMD = e and is called the epoch,* 

DMfi — nt^ 

TMA — a — longitude of the apse ; 
therefore, mean anomaly = AMfi — e - a, 
true anomaly = A MM' = TMM ' — TMA — 6 — a, 

12. 2h express the mean anomaly in terms of the true in 
a series ascending accordlny to the powers of e, as far as o’*. 

27r _ 27r • ^ 

periodic time h 

^ ^irh ^ h 

^irab 



* Wc shall sec (Art. 36) that the introductior of the epoch t is avoided 
in the Lunar Theory by a particular assumption ] but in tlie Planetary it forma 
one of the important elements of the orbit. 
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therefore h = wa* (1 — 6®)^, 

dt 

d6 h h * {1 -l-e cos (0 — aj}'* 

= ~ (1 — e^f {1 + e cos (0- a)}"® 

= ~ (I ~ \d^) {1 - 2e cos (^ - a) + 3e® cos® [Q — a)} 

= {1 — 26 cos (0 - a) + |e® cos2 (0 — a)} ; 

therefore + £ = 0 — 26 sin (0 — a) -f f 6® sin 2 (0 - a), 
or + £ — a) = (5 — a) — 26 sin (0 — a) + fe® sm2 (0 - a), 
the required relation. 

13. To express the true anomaly in terms of the mean to 
the same order of approximation. 

a=:/j^4.s- a-f 26 sin (0- a) — fe® sin 2 (0 — a)...(l) ; 

~ a = + £ - a first approximation. 

Substituting this in the first small terms of (1), we get 
0 — a = Mi + £ — a -h 26 sin [nt + £ — a)...a second approximation. 

Substitute the second approximation in that small term 
of (1) which is multiplied by 6, and the first approximation 
in that multiplied by 6®, the result will be correct to that 
term, and gives 

0 - a = Mi + £ — a -f 26 sin {wi -f £ - a + 2e sin (Mi + 6 — a)} 

- |e®sin2 (Mi + £ - a) 

= Mi + £ - a + 26 sin (Mi + £ — a) 

+ 46® cos(?^i^-£~-a) sin(Mi+£-“a)— |e® sIn2(Mi+£- a) 
= Mi 4-e — a4 26 sin(Mi + £ — a) + |e® sin 2 (Mi + £ - a), 
the required relation. 
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The development could be carried on by the same process 
to any power of e, the coefficient of e® would be found 

If sin 3 [nt + e - a) ~ sin {nt + £ - a), 
but in what follows we shall not require anything beyond 

Prol'lem of Three Bodies. 

14. In order to fix the position of the moon with respect 
to the centre of the earth, which, by means of the process 
described in Art. (9), is supposed brought to rest, we must 
have some determinate invariable plane passing through the 
earth’s centre to which the motion may be referred. 

If the sun and earth were the only bodies in the universe, 
then the plane Ox* the elliptic orbit (in which, according 
to the last section, the motion would take place) would be 
fixed, and might be taken as the plane of reference ; ^ but, 
as soon as wc take into account the action of the moon 
and planets — especially the moon — the plane ceases to be 
fixed, and some other plane must be found not affected by 
these disturbances. 

Theory teaches us that such a plane exists,^’ but as its 
exact determination can only be the work of time, the 


* See Pvinsot, “ Th'nAe ct cUUnnimtion de P rffdateur du systume soIah'Cy* 
where he proves that an invariable plane exists for the solar system, tliat is, a 
plane whose position relatively to the fixed stars will always be the same, 
whatever clianj^es the orbits of the planets may experience; but as its position 
depends on the moments of inertia of the sun, planets, and satellites, and ihei'e- 
fore on their internal conformation, it cannot be determined d jtriori, and ages 
must elapse before observation can furnish sufficient data for doing so d 
jwsterinri. 

This rcsidt Poinsot obtains on the supposition that the solar system is a 
fi’ee system ; but it is ijossible, as he furthermore remarks, nay })robable, that the 
stars exert some action upon it, it follows that this invnr-nhle plane may itself be 
variable, though the change must, according to our ideas of time and space, 
be indefinitely slow and small. 
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following theorem will supply us with a plane whose motion 
is extremely slow, and which may for a very long period, 
and to a degree of approximation far beyond that to which 
we shall carry our investigations, be considered as fixed and 
coinciding with its position at present. 

In what follows the action of the planets is too small 
to be taken into account. 

15. The centre of gravity of the earth and moon describes 
relatively to the sun an orbit very nearly in one plane and 
elliptic; the square of the ratio of the distances of the moon 
and sun from the earth being neglected.* 



Let 8^ JSj M be the centres of the sun, earth, and moon, 
O the centre of gravity of the last two. Now the motion 
ot O is the same as if the whole mass E-\-M were collected 
there and acted on by forces equal and parallel to the moving 
forces which act on E and M. The whole force on O is 
therefore in the plane 8EM\ join SG, 

Let L 80M— ct), and let m' be the sun’s absolute force. 

Moving force on in E8^ 

moving force on in M8, 


♦ This ratio is about and, as we shall see Art. (21), such a quantity 
we shall consider as of the 2nd^rder of small quantities, and its square therefore 
of the 4th order. Our investigations are carried to the 2nd order only. 
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These applied to G parallel to themselves are equi- 
valent to 


and 

but 


m'.E.QE m'.M.GM. 
— 


m'.E.SG ^m'.M.SG 


SE' 


SM^ 


MS-, 


E M _M+E 
GM~ GE ME ’ 


by the 
) triangle 
of forces. 


GM.GE 

whence E. GE= M. GM= (if+ E) . 

jxLjii 


Therefore substituting, and dividing by we get 

accelerating forces* 


m . 


and 

Again, 


OM*GE / I 1 \ . T. . riTijr 
— TiTn — 1 TTTta — direction GM. 

HE \SE^ 8W) ’ 

m\8G ( E ^ 

M^E \8E^^ 8MV 

8E^ ^8G^^rOE^^28G.GE cos a> 


0/72 iGE 

= 8G ( 1 4- — COSO? + 

\ bLr 


GE\ 

8 GV^ 


therefore 

similarly 


1 1 / 3GE 

JSE^ ~ SG^ V SG jj 

1 1 /, ZGM \j 

SM^ ~ SG^ V SG 


omitting quantities 
to be neglected. 


♦ In strictness it would be necessary, since we have brought S to rest to 
apply to both M and is", and therefore to G, accelerating forces equal and opposite 
to those which E and M themselves exert on IS ; but the mass of S is so large 
compared with those of E and that we may safely neglect these forces in 
this approximate determination of the path of G, the error being of a still 

higher order than that introduced by the neglect of 
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Tliereforc the accelerating force in the direction OM 


Sm\GM. GE 


8G^ 


cos (O 


== - 3 (accelerating force of sun on G) cosw 

= 0, according to standard of approximation adopted. 
And the accelerating force in G8 

m'.SG (E + M S{E.GE--M.GM) 

m . . , 

= to the same approximation. 


Hence, the force on is a central force tending to 8 
and varying inversely as the square of the distance ; there- 
fore the orbit of G about 8 is very approximately an ellipse 
with 8 in the focus ; and the plane of this ellipse is, as far 
as our investigations arc concerned, a fixed plane. 

This fixed plane is called the plane of the ecliptic^ or 
simply the ecliptic. 


16. We shall now suppose a plane through the earth’s 
centre parallel to the ecliptic; and this, when the earth is 
brought to rest, will give us the fixed plane of reference 
we require (14); the ecliptic then making small monthly 
oscillations from one side to the other of our fixed plane. 

17. Since G describes an ellipse relatively to the sun, 
the sun will describe the same ellipse relatively to G ; but, 
as seen from the earth, the orbit of the sun will be slightly 
different, as we shall now shew. 

First : The sun will have a latitude, that is, will be out 
of the plane of reference; and this latitude will be of the 
same name (north or south) as that of the moon, and 
deducible from it when the ratio of the distances of the two 




and moon are known. For if 8'EM' be the fixed plane 
through and 8\ G\ M' the projections of 8^ if, then 
8E8' is the sun’s latitude, and MEM' is the moon’s; and 
the two bodies arc obviously both on the same side of 


the plane. 

sin (sun’s lat.) = 


88^ _ 
ES^' 


GG' 

~E8 


sin (moon’s lat.) 


M EM . , , , , 

Now, it is known that ilf is about of E^ 

and i?lf of^^/Sf; 


therefore 


sin (sun’s lat.) = 


sin (moon’s lat.) 
—32400 


And as the moon’s latitude never exceeds 5° 9', the sun's 
latitude will always be less than 1". 

Secondly, as to the sun’s longitude: let ET be the 
direction of the first point of Aries, — that is, a fixed line 
in our plane of reference from which the longitudes of the 
bodies are reckoned. Then TE8' = 6' is the sun’s longitude 
as seen from E] but, as seen from G, the sun’s longitude 
would be the angle between G8 and a line through G (not 
drawn in the figure) parallel to ET. 

The difference between these two longitudes is the angle 

E8'G'] emd/if EG' 8' = (a, 

• -ncitm EG' . " sin'O . , , 

B,mES G =‘ 32 ^ approximately; 
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therefore sin ES'G' never exceeds ES^G' is a 

small angle not exceeding 7”. 

Also ES’ ^ S*G' < EG* < ott ^ • 

18. Now, by assuming the longitude and distance of the 
sun as seen from E to be the same as when seen from G^ 
we commit the above small errors ii the position of 8] 
that is, we assume the sun to be at 8'* instead of 8^ 8' 8’* 
being drawn equal and parallel to G'E, If our object were 
the determination of the sun’s position, it would be neces- 
sary to take this into account ; but when investigating the 
motion of the moon, the disturbing action of the sun will 
not, on account of his great distance, be appreciably altered 
by supposing him to be at 8” instead of 8 ; and the errors 
thus introduced are much too small to be considered. 

lienee we may assume that the motion of the sun about 
the earth at rest is an ellipse having the earth for its focus, 
and its equation 

u* = a {1 -t-e cos(0' — J*)}, 

where ^ is the longitude of the sun’s perigee and e = Is 
the eccentricity, and we are safe that no appreciable error 
will ensue in the determination of the moon’s place.* 


* That is, as far as the three bodies alone are concerned,- — but, since the 
attractions of the planets may, and iu fact do, disturb the elliptic orbit of 
the sun about 6?, the same cause will disturb the assumed orbit about E, A 
remarkable result of this disturbance is noticed in Appendix, Art. (108). 
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CHAPTER III. 


RIGOROUS DIFFERENTIAL EQUATIONS OF THE MOON’S MOTION 
AND APPROXIMATE EXPRESSIONS OF THE FORCES. 


19. The earth having been reduced to rest by the process 

described in Art. (9), let us take its 

centre E for origin. Let r, 6 be the 

polar coordinates of the projection M 

of the moon on the fixed plane of re- ^ 

ference, Q being the longitude TEM' measured from a fixed 

line Ex in that plane. Let s be the tangent of the moon’s 

latitude MEM\ so that MM' = rs. 

Next, let the accelerating forces which act Upon the moon 
be resolved into these three : 

P parallel to the projected radius ME and towards the earth, 
T parallel to the fixed plane, perpendicular to P, and i^ the 
direction of 9 increasing^ 

S perpendicular to the fixed plane and towards it. 

The equations of motion will be 


1 d 

r dt 

de ' 



20. These tlirce equations for determining the moon’s 
motion take the time t for independent variable ; but it will 
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be more convenient in the following process to consider the 
longitude as such, and our next step will be to change the 
independent variable from t to 6. 

From (i) we get 


,d9 

dt dt\ dt)~ dt ’ 


therefore 


Therefore 


= suppose, whence JBT ^ = 2V*. 

'~T. — 2 — ^ \t U — , 

at r ^ 

1 I f. ^ t T 




an important equation connecting the time t with the longi^ 
tude 0. 

. . dr dr dd „du 

^'■Z’de-dr-^di' 

d?r d I ^ydu\ o d 

df-d,{-^de) — ^‘'‘36V’Tii)- 

Subsatute in (ii), 

Iberefor, a-J, (;/*) + flV = P, 






and substituting for IP and for their values 2 J^d0 
T 

and ^ , then dividing by and transposing, 


d^ __ P T du 
AV hWdO 




this is called the diferential equation of the moon^s radius 
vector. 
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Lastly, 



from (4), F, - ifV + 

tberofore f>r - //V (g + «) + iff »• | ; 

therefore 


_Fs-8 


JT * 



this is the differential equation of the moon^s latitude* 


2^. If these three equations could be integrated under 
these general forms, the problem of the moon’s motion would 
be completely solved ; for as only four variables w, 0, 5, and 
t are involved (the accelerating forces P, P, and 8 being 
functions of these), the values of three of them, as w, 5, 
could be obtained in terms of the fourth t ; that is, the radius 
vector, longitude, and latitude would be known correspond- 
ing to a given time. 

But the integration has never yet been effected^ except 
for particular values of P, P, and 8 ; and the method which 
we are in consequence forced to adopt is that of successive 
approximation, by which the values of w, and 5 are 
obtained in series in which the terms proceed according to 
ascending powers of certain small fractions. Some one 
fraction being chosen as a standard with which all other 
are compared, the order of the approximation is esteemed 
by the highest power of the small fractions retained. 
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It Is usual to consider as a small fraction of the first order, 


consequently of second 

third 


and so on, other fractions being considered as of the 1st, 
2nd, &c. orders, according as they more nearly coincide with 
-ih-U} 

22. It is necessary therefore, before we can approximate 
at all, that we should have a previous knowledge (a rough 
one is suflScient) of the values of some of the quantities 
involved in our investigations; and for this knowledge we 
must have recourse to observation. 

We shall therefore assume as data the following results 
of observation : 

(1) The moon moves in longitude aoout thirteen times 
as fast as the sun. The ratio of the mean motions in longi- 
tude represented by m is therefore about and may be 
considered as of the 1st order.^ 

(2) The sun’s distance from the earth is about 400 times 
as great as the moon’s distance. 

Hence the ratio of the mean distances = is of the 
second order.f 

(3) The eccentricity e' of the elliptic orbit which the 
sun approximately describes about the earth is about 
and this, approaching nearer in value to ^ than to 

will be considered as of the 1st order. 

(4) During one revolution, the moon moves pretty ac- 
curately in a plane inclined to the plane of the ecliptic at 

* This approximate value of m is easily obtained j — the moon is found to 
perform the tour of tlie heavens, returning to the same position among the 
fixed stars, in about 27^ days; the sun takes 305 J days to accomplish the 
same journey. 

f The distances of the luminaries may be calculated from their horizontal 
parallaxes, found by observtt^ions made at remote geograpliical stations. (See 
the Author’s Astronomy^ chap. xvi.). 
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an angle whose tangent is about and therefore of the 
1st order.* 

(5) Its orbit in this plane is very nearly an ellipse having 
the centre of the earth in its focus, and whose eccentricity is 
about equal to our standard of small fractions of the 1st order, 
’"iz. ; and this will also be very nearly true of the projec- 

tion of the orbit on the plane of the ecliptic.f 

To calculate the values of the forces. 

23. Wc are now in possession of the data requisite for 
beginning our approximations, and we shall proceed to the 
determination of the values of T, and S in terms of the 
coordinates of the positions of the sun and moon. 

Let 8^ M be the centres of the sun, earth, and moon, 



E\ M\ the projections on the plane of the ecliptic, 
G the centre of gravity of E and M, 


* That the moon’s orbit during one revolution is very nearly a plane inclined 
as we have stated, will be found by noting her position day .after day among the 
fixed stars; and the sun’s path having previously been ascertained in a similar 
way, the rules of Spherical Trigonometry will easily enable us to verify these 
statements. {A stronomy, p. 92) . 

f The elliptic nature and the value of the eccentricity of the moon's orbit may 
be found by daily observation of her parallax, whence her distance from the earth’s 
centre may be determined : corresponding observations of her place in the heavens 
being taken, and corrected for parallax to reduce tliem to tlie earth's centre, will 
determine her angular motion. Lines proportional to the distances being then 
drawn from a point in the proper directions, the extremities mark out the form 
of the moon’s orbit. 

A similar method applied to observations of the diameter of the sun will 
determine the eccentricity of its orbit. {Astrouot.vt/, p. 138). 
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The forces we have to take into account are, the forces 
which act directly on if, and forces equal and opposite to 
those which act on E ] — these last being applied to the whole 
system so that E may he a fixed point. 


Attraction of E^ 


Forces on if. 
E 


ME^^ 


attraction of 8 
equivalent to ^ 


\n{Ma 

mr' * 


An ME 
.inifE 


W.SG , 
^j^3-pfirUo G8 


Forces on E, 

M 

Attraction of i/= , ,,^3 in EM 

ME^ 


attraction of S 

equivalent to 


vni %juj\jr , -riTiir 

'-77rSf- mEM 


\m^G 

83 


m'.8Q , „„ 

SE^ G8 


Therefore, the whole attraction upon M, when E is 
brought to rest, is 

E+M^ ,fMG EG\. ' ™ 

ME^ MV *" 

and n^8G ^^^5 ~ '8E^ parallel to G8. 


These are the rigorous expressions of the accelerating 
forces on AT, and they can be expressed in terms of the 
masses and coordinates of the bodies. 


Approximate values of P, T, 8, 

24. Our investigations will be carried only to the second 
order ; It will be sufficient therefore if, in the preceding, we 
neglect small quantities of the fourth and higher orders. 

Let /4-=P+ if, 

MGM' = tan"^5 = moon’s latitude, 

E'T the direction of the first point of Aries, 

8G = r' — = longitude of sun, 

M'E' = ^ = ^ ; FTE'M' = 0 = longitude of moon, 
/SP'JP = 0 - 0' = difference of longitude of sun and moon. 
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Now, SM^^SW+MM'^ 

= 8CP+ GM'^ - 28G. GM' cos 8GM’ + MW' 


therefore 


= 7» ; 


^QM' 


for 6’-~6* or 8E*M’ differs from 8GM* by less than 7', 
/ OM\^ . 

Art. (17), and f ^ is neglected, being of the fourth order. 


Similarly, 




therefore 


the accelerating forces on the moon are approximately 

+ ^5 ini direction ME^ 

and ~{GM'+ GE’) <ioa{0-0') parallel to G8', 

whence P= ME'^ coaMGM'- ^ M’E’ coa\0-0') 

= (1 - |s*) - ^ B + I cos2 [0 - 

u 

T= - ^ M'E' cos {0 - &) sin {0 - ff) 

_ Wl IM m t /% /W\ 

3C — ^ Sin 2 (0 — 0 ), 


S-{^. + %SlE)«nM(}M’ 

_ f /a m'r \/{l 4 s*)") s 

"Va +«•)■*■ /“ ; 


Vll + s'l 


O / ^ 9V ••• «> 

= /tw (s-|«) + -— ; 

Vi 


Ps-8 — ■ 


-B4-|cos2(0-^)}. 
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The differential equations in Art. (20), when these values 
of the forces are substituted in them, would contalu the 
variables u and but u is given In terms of 6' by the 
equation of Art. (18) u—a'{l+e cos(^' — f)}, and we shall 
find means to establish a connexion betVeen 6 and 6\ which 
will enable us to eliminate the latter. 


25. Before proceeding further, it will be important to 
consider the order of the disturhing 'effect of the sun’s action, 
compared with the direct action of the earth. If we examine 
the values of P, P, and >S, it will be found that the most 
important of the terms containing which are clearly the 
disturbing forces since they depend upon the sun, are in- 

7yiT 

volved in the form ^ , while those independent of the sun’s 
action enter in the form ^ . 

T 

We must therefore find the order of compared with ^ , 


or ot 

r 




Now the orbits being nearly circular, and m the ratio 
of the mean motions, Art. (22), we have 

mean motion of sun periodic time of moon 
.mean motion of moon periodic time of sun 

^ 2'7rr^ 27rr ^ 

therefore ^ ~ : 1 , 

T T 

or the disturbing force of the sun is of the second order 
when comnared with the direct action of the earth. 
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CHAPTER IV. 

DSTBGSATION OF THE DIFFERENTIAL EQUATIONS. 

Section I. 

Oeneral process described. 

26. If in the differential equations we substitute the 
values of P, 1\ and S just obtained, we are led to the 
following forms: 

which cannot be integrated ; but it is found possible to obtain 
a solution which will be a first approximation, if we retain 
only the important terms on the right-hand side of each 
equation. This first approximation being substituted in the 
equations, the terms next in order of importance may be 
retained and a new solution becomes possible, which in its 
turn is the basis of a further approximation, and so on. 

First, neglect the disturbing force of the sun which is 
of the second order, and also the moon’s latitude, which, 
as will be seen by referring to the expressions for the 
forces (Art. 24), will either enter to the Second power or else 
in combination with the disturbing force. 
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When this is done the equations become integrable, and 
values of u and 5 may be obtained in terms of 0 correct to 
the first order; this value of u will then enable us to get 
the connexion between 0 and t to the same order. 

[Let us, however, bear in mind that the equations thus integrated are not 
the diiferential equations of the moon’s motions, but only approximate forms 
of them ; and it is, therefore, possible that the results obtained may be ap- 
proximations to the true solutions only for a limited range of values. 

Whether they are so or not, will appear by comparing them with what we 
already know of the motion from observation ; and this previous knowledge^ 
in the event of their not being approximations, will probably suggest such 
modifications of them as will render them so.] 

To proceed to the next approximation: — We must sim- 
plify the right-hand members of our general equations by 
employing the values of u and s just obtained, retaining 
now terms of the second order which we had previously 
neglected. 

The integration can always be performed when the right- 
hand members of the equations are circular functions (sines 
or cosines) of 0 or of its multiples ; and as our first approxi- 
mation will give us the values of u and s in that form, the 
new expressions for the forces will also be of the same 
character, and the equations be again integrable. 

Thus we shall obtain new values of w, 5 , t correct to the 
second order. These values, introduced in the same manner 
in the second members and terms of the next higher order 
retained, will lead to a third approximation, and so on, to 
any order; except that if we wish to carry it on beyond 
the third, the approximate values of the forces, given in 
Art. (24), would no longer be suflSciently exact, and we 
must obtain more correct values from the expressions of 
Art. (23). 

27. We must here stop to notice a peculiarity in these 
equations, when solved by this process. We have said that 
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to obtain the values to any order, all terms up to that order 
must be retained in the second members ; but it may happen 
that a term of an order beyond that to which we are work- 
ing would, if retained, be so altered by the integration as 
to come within the proposed order. 

Such terms must therefore not be rejected, and we 
shall proceed to exa.nine by what means they may be 
recognised. 

28. Suppose then that after an approximation to a certain 
order, the substitutions for the next steps have brought the 
equation in u to the form 

+ + (?cos(2J^+^)+ , 

where the coefficient O is one order beyond that which wo 
intend to retain. The solution of this equation will be of 
the form 

w= -f G' cos(jp0+ jsr)+ 

Q' being a constant which may bo determined by substi- 
tuting this value of u in the differential equation, 

whence O' = ; 

1 -p ^ 

from which we learn that if p differs very little from 1, O' 
will be at least one order lower than O^ and will come 
within our proposed approximation ; and consequently the 
term O (io^[pd-^H) must be retained in the differential 
equation. 

The equation of the moon’s latitude being of the same 
form as that of the radius vector, the same remarks apply 
to it. 

29. If p is small, O' ^G very nearly, and the term 
G Qi}%{p6 -f- H) will not rise in importance in the value of w. 
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We should therefore be led to reject it In the differential 
equation, except for the following reason : — 

In finding the connexion between the longitude and the 
time (one of the principal objects of the Theory), we must 
use the equation 


dt _ 



Now, substituting for u its value in terms of 0, we shall 
find that such a term as O* Q,o^{p6-^n) in u gives rise to 
a term cos(^0 + //) in this equation, where G' and G' 
are of the same order, then 

G"co,[pd + E), 


whence t = — sin + 77 ) ; 

G*' 

therefore, when p is a small quantity of the first order, — 

will be one order lower than G'\ and the term will have 
risen in importance by the integration. 

T 

But yet further, if such terms occur in , they will be 
twice increased in value ; for they increase once in forming 
I j Sind once again, as above, in finding t. 


30. We have, therefore, the following rule: — 

In approximating to any given order to the values of u 
and 5, we niust^ in their differential equations^ retain periodical 
terms ONE ORDER leyond the proposed one^ when the coefficient 
of d in their argument^ is near unity. 


♦ The angle of a periodical term is its only variable part and is called 
’^he argument. 
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In approximating to the value of we must^ in the same 
equations^ also retain terms one order heyond the 2 )roposed 
one when the coefficient of the argument is near zero; and 

T 

when such terms occur in ^ they must be retained TWO 

h u ^ 

ORDERS beyond the proposed approximation.^ 


31. We shall here, for convenience of reference, bring 
together the equations and the approximate expressions for 
the forces : — 


^ ^ T du ^/(Fu \ [T , ,, 

/tV ^ Jav 


(Tu 


d& 


d's . Ps-8 


dtr 

dt 


• dd 
T 


P 

ikJ 


Ps-S 


|-a(‘+^/a- 4 ‘ 

% + 3 cos2 (0 - 6 ')} (2?), 




3 TwV® . ^ 


?5^1l+cos2(9-e')) IF). 


* Instead of the forces which really act on the moon, we originally substituted 
three equivalent ones, P, Tj S j these again are, by the preceding expressions, 

p 

replaced by a set of others. For, we may conceive each of the terms in , <fec., 

to correspond to a force, — a component of P, Ty or 8; each force having the 
same argument as the term to which it corresponds, and therefore going through 
its cycle of values in the same time. Now, by Art. (29), when the coefficient of 
6 in the argument is near unity, the term becomes important in the radius vector, 
and when near zero, in the longitude : hence, a force whose period is nearly the 
same as that of the moon, produces important effects in the radius vector ; and a 
force whose period is very long will be important in its effect on the longitude. 
See Airy's Tracts, Planetary Theoryy p. 78. 
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Section 1 1. 

To solve the Equations to the first order. 

32. Neglect all terms depending on the disturbing force, 
f.e. those which contain m] such teriiis being of the second 
order. Art. (25). 

The latitude s of the moon can n^'ver exceed the inclina- 
tion of the orbit to the ecliptic ; but this inclination is of the 
first order, Art. (22), therefore s is at least of the first order 
and s* may be neglected. 

WL P T ^ Ps^S ^ 

Whence ~ = 0; =0, 

and the differential equations become 
d^u fi 


^ n 


whence {1 + e coa(0 — a)} ; or, writing a for ^ , 


u = a {1 cos(0 — a)} (JJJ, 

and 5 = A sin (0 — 7) ( 8 ^)y 


€j a, kf 7 being the four constants introduced by integration. 

33. These results are in perfect agreement with what 
rouffh observations had already taught us concerning the 
moon’s motion Art. (22); for w = a {1 +e cos{0— a)} repre- 
sents motion in an ellipse about the earth as focus. 

And s = A; sin (0 — 7) indicates motion in a plane inclined 
to the ecliptic at an angle tan'”^Ar. 

For, if TOM' be the ecliptic, 

M the moon’s place, MM' an arc 
perpendicular to the ecliptic, then 




SOLUTION TO FIKST OliDCK. X) 

TM'^6] and if TO be taken equal to 7, and Oil/ joined 
by an arc of great circle, we have 

tanil/Jf ' = taniliOd/' sin OM ' ; 
or s = tan M OM ' sin ( ^ - 7), 

which, compared with the equation above, shews that 
ilfOii' = tan“U*. 

Therefore the moon moves in a plane passing through 
a fixed point 0 and making a constant angle with the 
ecliptic. 

The values of a and 7 introduced in the above solutions 
are respectively the longtitude of the apse and of the node. 

34. What the equations cannot teach us, however, is 
the magnitude of the quantities e and k. For this we must 
have recourse to observation, and by referring to Art. (22), 
we see that e is about and k about that is, both 
quantities are of the first order. Their exact values cannot 
yet be obtained: the means of doing so from multiplied 
observations will be indicated further on. 

35. Lastly, to find the connexion between t and 0, the 
equation ((7) becomes, T=0, 

^ 1 _ 1 1 

dtf hu^ lia^ [1 + e cos(t/ ~ a;}* ‘ 

Now this is the very same equation that wo Lad found 
connecting t and 0 in the problem of two bodies^ Art. (12), as 
we ought to expect, since we have neglected the sun’s action. 
Therefore, if p be the moon’s mean angular velocity, we 
should, following the same process as in the article referred 
to, arrive at the result 

0—pt + £ + 26 8in(^^ + £ - «)+ fe* sin2 + £ — a) + , 

which is correct only to the firs* order, since we have 

F 
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rejected some terms of the second order by neglecting the 
disturbing force. 

36, The arbitrary constant e, introduced in the process 
of integration, can be got rid of by a proper assumption: 
this assumption is, that the time t is reckoned from the 
instant when the 7nean value of 0 is ^jcro.* 

For, since the mean value of 0, found by rejecting the 
periodical terms, is + if, when this vanishes, ^ = 0, we 
must have e = 0 ; therefore 

d-pt-\‘2e sin(jp^ — a) 0^, 

correct to the first order. f 

37. We have now obtained three results, 8^, 0^, 
as solutions to the first order of our difftrential equations, 
and we must employ them to obtain the next approximate 
solutions; but, before they can be so employed, they will 
require to be slightly modified — in such a manner, howevsr, 
as not to Interfere with their degree of approximation. 

The necessity for such a modification will appear from 
the following considerations : 


♦ When a function of a variable contains periodical terms which go through 
all their changes positive and negative as the variable increases continuously, 
the 7/iean value of the function is the part which is inde|)endent of the periodical 
terms. 

t We shall also employ thi.s method of correcting the integral in our next 
approximation to the value of tl in terms of t ; and if we purposed to carry our 
approximations to a higher order than the second, w'e should still adopt the 
same value, that is, zero, for the arbitrary constant introduced by the integration. 
To shew the advantage of thus correcting with respect to mean values : suppose 
we reckoned the time from some dejinite value of 0, for instance when 0 = 0; 
then, in the first approximation, 

0 = 1 + 2e sin (• — a) 

is the equation for determining the constant i, and in the second approximation, 
c would be found from 

0 = € + 2c sin (« - o) + sin 2 (• -a) + , 

giving difiterent valaes of « at each successive approximation. 
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Suppose we proceed with the values already obtained; 
we have, by Art. (24), 

P 


AV 2AV *“•' 


= a(l - (1 + e cos(^-o)}'’ 


= a + -4 co8(0 - a) + ; 

and this being substituted in the differential equation {A) 
of Art (31 j, gives 

^ + M = a4- + ^ cos(^-a)+ , 

the solution of which is 

tts=:a {1 +ecosC0 — a)} 4* + sinf^ - a). 

Our first approximate value ?^ = a {1 + c co8(^ - a)] is here 
corrected by a term which, on account of the factor admits 
of indefinite increase, and thus becomes ultimately more 
important than that with which we started. Such a cor- 
rection is inadmissible, for our first value is confirmed by 
observation (22J to be very nearly the true one. The moon’s 
distance, as determined by her parallax, is never much less 
than GO times the earth’s radius; whereas this new value 
of when 6 is very great, would make the distance in- 
definitely small. On the same principle, we sec that any 
solution, which comprises a term of the form sin(0 — a), 
cannot be an approximate solution except for a small range 
of values of 0, 


Such terms, ‘ if they really had an existence in our system, must end in ita 
‘ destruction, or at least in the total subv^ersion of its present state ; but when 
‘they do occur, they have their ongin, not in the nature of the differential 
‘ equations, but in the imperfection of our analysis, and in the inadequate repre- 
^sentation of the perturbations, and are to be got rid of, or rather included in 
* more general expressions of a periodical nature, by a more refined investigation 
(than that which led us to them. The nature of tl is difficulty will Ijc easily 
‘understood from the following reasoning. Suppose that a term, such as 
‘a sin (.id 4 - 3), should exist in the value of a, in which A being extremely 
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‘ minute, the period of the inequality denoted by’ it would be of great length ; 

‘ then, whatever might be the value of the coefficient a, the inequality would 
‘ still be always confined within certain limits, and after many ages would return 
< to its former state. 

*Supi)oso now that our peculiar mode of arriving at the value of u led us 
‘to this term, not in its real analytical form a sin ^ 10 + /^), but by the way of 

* iU development in i)Owers of 0, a -f /i0 + X0'^ + A'c. ; and that, not at once, but 
‘ piecemeal, as it were ; a first approximation giving us only the term a, a 
‘second adding the term /ffi, and soon. If we strapped here, it is obvious tha^ 
‘we should mistake the nature of this inequality, and that a really periodical 
‘function, from the effect of an imperfect approximation, would appear under 

‘ the form of one not periodical These terms in the value of t/, when they 

‘occur, are not superfluous; they are e.sseiitial to its expre.ssion, but they lead 
‘us to ciToneous conclusions as to the stability of our system and the general 
‘ laws of its perturbations, unless we keep in view that thnj are onhj jmrts of 
^ series \ the principal parts, it is true, when we confine ourselves to intervals 
‘ of moderate length, but which cease to be so after the lapse of very long times 

* the rest of the series acquiring ultimately the jireponderance, and compensating 
‘the want of periodicity of its fiist terms.-Siu John Herschel, 
AhiropoUtana—ViiY Astuonomy, p. 079. 


38. To cxtritcate ourselves from this difficulty, and 
to alter the solution so that none hut periodical terms 
may be introduced, let us again observe that the equation 

= ^ = which gave the solution U and thus led 
du h 


to the difficulty, is only an approximate form of the first 
order of the exact equation. Any value of w, therefore, 

which satisfies the approximate equation + ^ = « fbe 


first order, and which evades the diflSculty mentioned above, 
may be taken as a solution to the same order of the exact 
equation. 

These conditions will be satisfied if we assume 


a = a {1 + e cos(c0 - a)} ; 
then 4- M = a + 0^(1 - c*) cos(c0 — a) 

= a to the*first order, 

provided 1 — c* be of the first order at least. 
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39. The Introduction of the factor c in the argument is 
an artifice due to Clairaut; and Its effect, as we shall see 
in Art (66) is equivalent to supposing that the apse of the 
moon’s orbit Is not fixed. Now, although the observations 
recorded in Art. (22) did not suggest this, we must bear in 
mind that they were extremely rough, and carried on only 
for a short interval; but when they are made with a little 
more accuracy, and extended over several revolutions of the 
moon, it is soon found that both her apse and the plane 
of her orbit are in constant motion. 

Clairaut was fully aware of these motions, and there is 
no doubt that he was led to the above form of the value 
of u by that consideration, and by his acquaintance with the 
results of Newton’s ninth section, which, when translated into 
analytical language, lead at once to the same form.** 

We might, therefore, taking for granted the results of 
observation, have commenced our approximation at this step, 
and have at once written down u = a[l+e co8(c0--a)} ; but 
we should, in so doing, have merely postponed the difficulty 
to the next step, since there again, as we shall find, the 
differential equation is of the form 


dO'^ 


■f u = a function of 


♦ Kewtou lias there shewn, that if the angular velocity of the orbit bo to 

O’t » 

that of the body as G — F to G, the additional centripetal force is — 
the original force being Therefore 




-f « = 


M 6^2- F2 


dhi F"^ _ fx F^ fi _ 


F. 


« = a |l 4* c cos 0 — , 


where is the same as our c, 

6r 
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the correct integral of which would be 

M = .4coa(0- jB) + , 

and this would at the next operation bring in a term with 
6 for a coefficient, which we now know must not be. We 
ehall, therefore, hereafter omit such terms as A cos(0 — -Bj 
altogether, and merely write 

ie = a {l + eco3(o0~ aj} -f 

40. The value of c will be found more and more' cor- 
rectly at each successive approximation, by always assuming 
these terms as the first two terms of the value of w, then 
substituting in the ditferential equation and equating coeffi- 
cients. It will thus be found that ae(l — c*) must equal the 
coefficient of cos(c^ — a) in the differential equation; and this 
will enable us to determine c to the same degree of approxi- 
mation as that of the differential equation itself. See 
Arts. (52) and (103). So far, all that we know about c is that 
it differs from unity at most by a quantity of the first order. 

41. In carrying on the solution of s, the same difficulty 
arises as in w, and it will be found necessary to change it into 

s — k — 7), 

g being a quantity which differs from unity at most by a 
quantity of the first order. 

The introduction of g is connected with the motion of the 
node in the same way as that of c is with the mqtion of the 
apse (85j ; and the value of g will be determined by a process 
similar to that explained above for c. Arts. (51) and (102). 

42. The connexion between 6 and t will also be modified 
by this change in the value of m, 

dt _ 1 1 

dQ ^ Aa* {1 + e cos(ct/ — a)]’* ’ 

d.ci _ J 1 

d.cQ ” Aa* {I cos(c^- a)}* * 


or 
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Here c0 and ct hold the places which 0 and t occupied 
in (35); therefore 

c0 = cpt + 2e sin(c/?^ — a), 
or 0=pt + 2e sin (cpt — a), 

• e 

+0 the first order, since - = e to the first order. 

^ c 

43. Since the disturbing forces are to be taken Into 
account in the next approximation, we shall have to- use 
the value of u found in (18), which is 

= {l+e'cos(fl'-?)}: 

but this introduces ff ; we must therefore further modify it 
by substituting for 0* its value in terms of 0^ and it will 
be found sufficient, for the purpose of the present work, 
to obtain the connexion between them to the first order, 
which may be done as follows : 

Let m be the ratio of the mean motions of the sun and 
moon, 

p\ p their mean angular velocities; .\p —mp^ 

P^pt mean longitudes at time /3 being the su’^’s- 

longitude when ^ = 0, 

0\ 0 true longitudes at time f, 

a longitude of perigees when 0;. 

therefore O' - ?= sun’s true anomaly, 

and - mean anomaly. 

But, by Art. (13), 

true anomaly = mean anomaly + 2e sin (mean anomaly) + &c. ; 
therefore O' = p't-^ ^ 4- 2e' sin ( + /3 — f ) + 

= mpt-\ ^ +2e sin + •••••• 

= + +2e'sin(m0 -f yS- f) 

to the first order ; 
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because pt^d — 2e sin (c0 — a) to the first order by (42). 
Whence u = a [I + e cos (jn0 + — f)} to the first order. 

44. The values of sin 2 (0- 0') and cos 2 (0-0') can also 
be readily obtained to the same order ; 

sin 2 (0-~ff) 

= sin {(2 — 2m) 0-2/3 — 4e' sin [m0 + ^ — f)} 

= sin{(2 — 2 >/ 2 ) 0 — 20} —ie s\n{m0 /3 — f) ^o^[{2-2m)0 —20] 
= sin {(2 - 2m) 0 — 20] — 2e sin {(2 — m) 0-0 - 

^2e sin [(2 - 3m) - 30 + ^}. 
Similarly, cos 2 (5 — 0') = cos {(2 - 2?w) 0 — 20] 

— 2e' cos{(2 — m)^- +2e' cos {(2 ~ 3m) 0 - 30 -i- ^]. 

The first term of each of these is all we shall require. 

Section III. 

2h solve the equations to the Second Order, 

45. Let us recapitulate the results of the last approxi- 
mation. 

u — a {1 + e cos(c0 - a)}, 
u = a' {1 + e' cos(w0 0 - f )|, 
s — Jc sm(g0 — 

0 — ff — (l — m)0-0—2e 8in(m0 + iS— f). 

These values must now be substituted in the expressions 
for the forces, retaining terras above the second order, when, 
according to the criterion of Art. (30), they promise to be- 
come of the second ord^^r after integrating. 
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The differential equations (^) and (B) will then assume 
the forms 




'd&‘ 

and the integration of these will enable us to obtain u and $ 
to the second order; after which, equation [C) will give the 
connexion between 6 and t to the same order. 


46. The quantity which we shall meet with as a 

coefScient of the terms due to the disturbing force, can bo 
replaced by m being the ratio of the mean motions of 
the sun and moon. 

So long as we neglected the disturbing force, h and a 
had determinate values; — they belonged to the ellipse which 
formed our first imperfect solution, and would therefore bo 
known from the circumstances of motion in that ellipse at 
any instant; h being double the area described in a unit 
of time, and a the reciprocal of the semi-latus rectum. It 
would consequently be impossible to assume any arbitrary 
connexion between them. But, when we proceed to a second 
approximation and introduce the disturbing force, there is 
no longer a determinate ellipse to which the It and a apply : 
the equation y, = K^a of Art. (32) merely shews that a and h 
must refer to some one of the instantaneous ellipses which 
the moon could describe about the earth if the disturbance 
were to cease, and we are at liberty to select any one of 
these which will allow us to proceed with our approximation. 
The particular ellipse is determined by the above assumed 

TYiCl^ 

relation and the selection is suggested and 
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justified by the following reasoning: 

•p average period of moon about earth 
p period of sun about earth ’ 

but, since the instantaneous ellipses are nearly circles, we 
have, as in Art. (25), 

(period of moon about earth)^ _ md^ . 

(period of sun about earth)* pd ^ ’ 

therefore if a be properly chosen, 

2 _ md^ _ rdd^ 

^ pd dd * 

47. If we examine the equations of Art. (31) it will be 
seen that we cannot solve the differentia^ equation in n to 
the second order until we know the value of s accurately to 
that order ; for, suppose 

s A sinC^0--7)-f Fj, 

where contains the terms of the second order not yet 
found; and suppose that among the terms in F^ there be 
on.> iisin(Z0 — X), then this value of s substituted in the 
expression {B) gives 

~ ^ sin [<7^- 7) + L sin(Z0 — X)...]*} — ... 

= a U - sin*(;70-- 7) - ZlcL sin(p0 — 7) sin (75 — X)...}, 

and equation (A) becomes, so far as this term is concerned, 

i + «) = 1 +. ..- ^kL [cos ((y - Z) ^ - (7 - X)}J 

+ fA-i [cos {f^r + Z) ^ _ (7 4. X) j]. . . . 

Now, if it should happen that Z is a small quantity, then 
the two coeiBcients g-.l and ,7 + Z will be both near unity. 
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and therefore these two terms, although of the third order 
in the differential equation, would become Important in the 
value of u. Again, if I is near 2, then l-g will be near 
unity, and the same remark will apply. 

Such terras do not really occur in the value of 5 , but 
there is none the less a necessity for verifying the fact. 

If I is near unity, then the coefficient g — l becomes small ; 
and although the corresponding term in u is still of the third 
order, it becomes of the second order, and therefore impor- 
tant in Art. (30). 

If I is not near 0, 1 or 2, the term will not affect, to 
that order, the values of u or of t. 


48. Again, it is necessary to have the value of u correct 

T 

to the second order, before we can compute to the fourth 

order, with a full assurance that we have secured all terms 
of that order in which the coefficient of the argument 
is nearly zero*; — such terms becoming, as we have seen 
Art. (30), important in the value of t. 

For, assume that 


u = a {1 + c cos[c9 — a) + TFJ, 



3nia^ {I +e cosfm0+/8- f)}* 
2AV {1 + e cos {cO - a) + 


sin {(2 — 2m) 0-2^} 


= —|m’*(l...— 4irj sin {(2- 2m) 0-2/3}. 

Therefore if contain a term It cos [r0 — p) of the second 
order where r is near 2 (except 2 - 2m) there will be a 


♦ This was pointed out by Mr. Walton in the Quarterly Journal of MathematicSf 
Vol. IX., p. 227 ; and in this edition the text had been modified in accordance 
with his suggestion j but the method given by him fo'* finding the value of u is 
defective on account of its not mentioning the necessity for the previous deter- 
mination of 8 . 
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T 

corresponding term in ^3 of the fourth order with a small 

coefBcient of and therefore an important term in the 
value of t. 

49. For the sake of simplifying the expressions we shall 
generally omit a, and 7 in the argi.ments. If we remark 
that c and a always appear together in the form — a, it 
will be sufficient to write c6^ Instead of cd — OL] and, for 
the same reason, we may replace m6 + ^ and gO — y by 
m0^ and The suffix Q will remind us of the omission, 
and we may at any stage re-introduce the omitted symbols. 


To compute s to the second order. 

50. Since we must first determine the value of 5 , it 
becomes necessary to examine the differential equation [B) 
to see whether any difficulty of a similar kind to those we 
have been here considering may arise, or whether the 
elements on which the computation depends are already 
known with sufficient accuracy. We fortunately find that 
all the terms of the right-hand member of the equation are 
of the third order, and therefore that the values already 
obtained will ensure accuracy to the second order, and we 
may safely begin our work here. 


51. To the first order we have 
s=^k sin^0„ 
ds 


(because 5 ^ = 1 


= first order, 
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. Fs- 8 

In the expression tor - ^—3 — , which Is of the third order, 
it will be suiSScient to put a for u\ and a for w, then 

= - 2 (1 + cos (2 - 2m) 0,} 

= - sin (2 - 2 m --g) 9 ^^ 

the other term of the third order being rejected because the 

coefficient of 6 Is not near unity. 

mi p • T , ds , 

ihe next term or equation (i>j 

T 

of the first order, it will be suflScient to have correct 
' h u 

to the second order, 


AV ''de^~ 

= — \w?k sin (2 — 2m — g) 9 ^^ 

neglecting the other term of the third order because the 
coefficient of 0 is not near unity. 

fd^s \ [ T 

The last term of equation {B) is 2 f ^ 

r T 

which =0 to the third order, since h^jd 9 depending on 
the disturbing force is of the second order, and the other 

To solve it, assume 
h sinpi?^ 

-f A sin (2 - 2m - 5^) 9 ^. 


factor = 0 to the first order. 
The equation in $ becomes 


d^s 


— \yr^k sin5^6^j 
+ |m’A sin (2 - 2m - 9 . 
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Substitute, and equate coefficient of like terms, 

A (1 _ - ^m‘k ; jr = (1 + = 1 + f m’, 

A[l-{2-2m-g)‘] = %m‘k ; A = fZ 

Therefore, accurately, to the second order,* 

« = &sm(^0 — 7) 4 ^mk sin {(2 — 2m ~ (^r) 0 — 2^9 + 7}.. 

The terra of the second order, having a coefficient near 
unity, need not be taken into account in finding the value 
of u (47) ; but it will have to be brought in subsequently 
to determine t 


To determine u to ike second order, 

52. We must compute the right-hand member of the 
equation (A), Art. (31), accurately to the second order, and 
include those terms of the third order which have a coeflScIent 
near unity. 

To the first order we have 


tt = a(l + ecosc^j), 


^ = — sin c0^ = — ae sin cd^ 


because c = 1 
to first order, 


d^u 

jp+u-a. 


From the expressions {D) and (-ff), Art. (31), we get 


♦ No complementary term P cos (0~Q) is added; for, though by the theory 
of differential equation? this would form a necessary part of the solution, we 
have seen Art. (39) that it cannot, in this shape, form a part of the correct value, 
but wiU be comprised in the terms whose argument is ^0 y. 
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- {1 + 3e cos(w20j — f) — 3c cosc0J {14-3 cos(2 — 2m) 

a = 1 — 1^* 4- f /iJ* cos cos (2 — 27?i) 

4- f wi*c coso^j 4* '4^*6 cos (2 — 2?n - c) 6^. 
The other terms of the third order, which arise in the 
development of the expression, are neglected because the 
coefficients of 6 are not near unity, 

^ = _ 3mV^ {l + f'co8( >ng. - ^ 

h'u^ 2/iV(l + ecosc^,)* ' ^ * 

= - [1 4- 3e cos [m6^ — f ) — 4c cosc^J sin (2 - 2m) 6^ 

= — sin (2 — 2m) 6^ 4- 3m*c sin (2 — 2m — c) 

AV ' S “ 


= f ?/i*ae cos (2 ~ 2m — c) 0, ; 

rejecting the other terms of the third order according to (30). 

f iTs (2 ~ 2m) ^ cos (2 — 2m - c)^ 

Jhu ^2- 2m ^ ^ * 2-2m-c ^ ^ ‘ 

= I cos (2 - 2m) - 3m*e cos (2 — 2m - c) 

m*a cos(2-2m) 0^-6m*aeco8(2-2m~c)0,. 


The equation in u becomes 


1 /ePu 


a \d& 


I* “*■ 


„) = 


1 - — 

4* Im'^c COSC0J 

4 - f A* COS 290 J 

— 3m* cos (2 — 2m) 0^ 

4* ym*e cos (2 - 2m - c) 0^. 


To integrate this, assume 
u _ 
a 

l+A 

4-ecosc0, 

4-460082^0, 

4- (7cos( 2 — 2m) 0, 

4- D co8(2 - 2m — c) 0,. 
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Substitute, and equate the coefficients of like terms : 


c(l — c") = 


C {1 - (2 - 2>n)’} = - 3m’, 


Z) {1 — (2 — 2m — c)’} = \fw?e. 
Therefore 




7?=2.i.;^=-jA’, 


(7= 


3>k’ 

l-'i 


+ m’, 


Z)= =Uwe. 


-f- e cos(c^ — a) 
u — ai — i/t:^cos2 — 7) 

+ cos ((2 ~ 2/??) 0 - 2^} 

, + V 1(2 — 2m — c) 0 — 2/3 + a} 
This Is the accurate value of u to the second, order. 


U. 


To find t to tlui second order. 


53. The equation to be considered is 


dj 

dd 



and, making use of the values of s and u obtained to the 
second order, we must retrace our steps to recompute u and 
T 

yj -g so as to include those terms of the third order in u in 

\h 

•which the coefficient of the argument is near zero, and 

T 

similar terms to the fourth order in Art. (48). To do 

hw ^ ‘ 

this we need not repeat the whole work, but only so much 
of it as will enable us to pick out the additional terms in 
question. 
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1 P 
a Vu 


~ = 1 — ^ {^• sin 4- ^mh sin (2 — 2 w 2 — g) 6^ 


— |?72’{l4-3r/ co^(m0-^)—3c cos^?^J{l+3 gos(2-2/2/)0j} 
= 1...— grf sin< 7 ^^ sin (2-2?>i-//) 0- coii(7n0^—^) 

= 1...— -ru'^nk^ cos (2 — 2?« — 2g) 0^ — cos (7n0^ - f ), 
T ^ 3m (I + <;' cos??i^j - sin (2 — 0^ 

2lca^ ( i + cos c0^ — \/v^ cos 2(/0^ )* 

= — I sin (2 — 27n) [I + e cosc^^ - :{/,•’* cos 2 ,( 7 ^j} 

= - Im^ siu(2 — 22 ?i) 0^ {(l + c cosc^^)"‘ + /•:* cos2(y^J 
= — sin (2 — 2m) 0^ {...10c* cos‘^c0j + /j* C082^^j} 

= — |w 2 * sin (2 — 2)n) 0 ^ {...5c* cos2o0j + /:* cos 2^0, j 
= — sin(2-2/?i-26*)^^~ sinf2 — 2 ?? 2 ~ 2 ^) 0 ^. 

These are the only important addillonal terms. 


r T 

hi^u^ (2— 2m -2c) 0^— gm/u* cos (2- 2m -2^) 0^^ 


because 2 -2?n - 2c = - 2??2 + and 2 - 2?n -2g = - 2m — 
Substituting in the dilFcrcntial equation for w, the neto 

terms in ^ will be of the fourth order and must bo 

th u dt/ 


neglected, but 
terms 


2 



d0 will bring in the two 


— cos (2 — 2m — 2c) 0^ and -> ^m]da cos (2 — 2m — 2g) 0 ^ ; 

and we get, so far as important additional terms are concerned, 


dy 

d&^ 


H- ttj = 1...+ '^rne^ cos (2 — 2m — 2c) 0^ 


4- (~ f% 4- 1) mF cos (2 - 27?i - 2g) 0^— |mV cos (m0^— f), 
and the corresponding terms in u will also be of the third 
order, viz. 

~ =. . . ^7ne^ cos (2 - 27w — 2c) 0^ 4* Y\mA;* cos (2 - 2m - 2r/) 0^ 

— |7>iV cos(?n0j — S'). 


II 
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Tho complete value of u therefore, so far as Important terms 
are concerned for determining the value of is 

- = 1 - + e — JF coB2gd^ 

ct 

+ m* cos(2 — 2m) 0^ + cos(2 — 2m — c) 0^ 

H- cos (2 - 2m - 2c) 0^ + cos (2 - 2m — 2g) 0^ 

— f mV cos(m0j, — f), 

r T 

and the complete value of d0 will be 
r T 

J^g d0=^m^ cos (2 - 2m) 0^ — ^/m* cos (2 - 2m — 2c) 0^ 

— §mi** cos (2 — 2m — 2^) 0^^ 

the other term of the third order has no influence on tj the 
coefficient of 0 not being small. 


54. We are now in a position to find the value of t to 
the second order : 


d0 ■" M 




If- =1 + ^00300, + W where W contains all terms of the 
a ’ 

second and third orders in the expression above, then 

= ^, {1 — 2 (e cosc0,+TF) + 3(e cosc^j+TT)*— ...} 

= ^, (1 — 2e cosc^j — 2W+3e‘ cos‘c0, + 6eW cosc^,)... 

= (1 + |e’ - 2e 5OSC0,- 2 TTh- f e’ cos2c^,+ 6e W coacO^). 
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The only part of ^eW cosc0, which gives a term of the third 
order in which the coefficient of the argument is small, is 
COSC0J {...+ ^me cos(2-- 2m — c) 0^ 
of which one term is 

+ ^^md* cos(2 — 2m — 2c) 6^. 

Therefore ^ is equal to 

' 1 + |c* + f A'* + — 2c COSC0, + f c* cos2c9^ 

1 + i A* cos 2g9^ — 2m* cos (2 — 2m) 0^ 

liG? — cos (2 — 2m - c) 0, -f 3m V cos {m9^ — f ) 

^ me* cos (2 — 2m — 2c) 9^ cos (2 — 2m — 2^) 0^. 
Also from above 

fT |m* cos (2— 2m) + ^mc* cos (2 — 2m— 2c) 

j + §mi* cos (2 — 2m — 2g) 9^. 

We have now to multiply these results together, and we see 
that the terms having for coefficients ^mc* and §m7c* will 
disappear in the product.* These were the terms which 

T 

were originally of the fourth order in . 


14 |c* 4- |A* 4 m* — 2c cosc0, + |c* co32c0j 
4 \k^ cos2^0j - \^m* cos (2 — 2m) 9^ 

— ^mc cos (2 — 2m — c) 4 3m*e' cos(m0j — ?). 


* The principle of the superposition of small disturbances might have led us 
^ priori to anticipate the disappearance of these terms : — ^for, when approximating 
to the second order we may, according to this principle, calculate separately the 
effect of the disturbing action of the sun, which is of the second order, and that 
of any other independeni disturbance of the same or of a higher order, and then 
add the results. No terra therefore which bears a trace of the action of both 
these forces should present itself to the second order. 

Now since neither k nor e enters into the argument of any term, it is obvious 
that a terra whose coefficient contains or e* at any stage will retain thena 
through all subsequent operations ; and if the terras having (2 -- 2wi — 2^) 
or (2 — 2 /» — 2 tf) 01 for argument could remain in the result to the second order 
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Let 






therefore p = Tic? (1 ~ ^e* — |7c* — n?') to the third order ; 
therefore, multiplying by p and integrating, we get, still to 
the second order, 

pt=0—2e sin(c0— a) + S6* sin2 (c9 — ol) 

+ sin 2 (/7^ — 7) 

«- ym^sin{(2 -27w) 0 — 2/3} 

- y m<3sln{(2— 2 w2-c) 0-2y84-a} 

+ 3 W sin (jnO + /3 — f ) 
no constant is added, the time being reckoned from the 
instant when the mean value of 6 vanishes, for the reasons 
explained in (Art. 36). 


Coordinates expressed in terms of the time. 

55. The preceding equations 8^^ 0.^ give the reci- 
procal of the radius vector, the latitude and the time in 
terms of the true longitude; but the principal object of the 
analytical investigations of the Lunar Theory being the 
formation of tables which give the coordinates of the moon 
at stated times, wc must express le, 5, and 6 in terms of t 

To do this, we must reverse the series pt=d — &c., and 
then substitute the value of 0 in the expressions for u and $. 

Now 0=pt + 2e sin (c0 - a) to the first order 
sin (cpt — a) ; 


they would there have /i;* or ^ rcsijectively in the coefficient, since they have 
them when they first appear. But i-, and m being independent, such terms 
must have arisen from tlic combination of others wliieli originally (in the 
expressions of the forces) involved these quantities separately: — terms like 
cos2(^t) — y) and e-cos2(c0 — a) combined with terms involving m: that 
is, terms which are representatives of disturbing forces of the second order (see 
note p. 31) combined wi^'h others depending on the sun’s action. Such terms 
wust therefore be of an order beyond the second. 
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therefore c0 — a = — a + 2e sin {c.pt ~ a) to the first order, 

2e sin (cd — a) = 2e {sin (cpt — a) + 2e sin (cpt — a) cos (cpt ~ a)} 

to the second order? 
= 2(3 sir (cpt — a) + 2e® sin 2 (cpt - a) ; 

and as 0 and pt differ by a quantity of the first order, they 
may be used indiscrimirately in terms of the second order; 
therefero 

sin(cp^ — a)+ sin2 (cpt — a) k 

— sin 2 (gpt — 7) 

+ y sin {(2 — 2m) pt - 2/3} . . . 

+ V {(‘^ 

— sin (mpt + /8 ~ f) / 

56. In the value of u given in Art. (52), substitute pt for 
6 in terms of the second order, and pt + 26 sln(c^^ - a) in the 
term of the first order ; then 

1— 6*+6 cos(c^^— a) +6® cos2(cp<— a) 

— cos 2 (gpt - 7} 

+ cos {(2 — 2m) pt — 2/3} 

4- cos{(2— 27/i-c)2^^—2^-4-a} 

57. Similarly, the expression for s becomes 
s = A sin [(gpt — 7) 4 26 sin (cpt - a)} 

4- %mk sin {(2 — 2m — ,<7) - 2/3 + 7} ; 

^m(gpt- 7) 

+ 6 sin {(5^ •\-c)pt—0L — 7} 

- 6 sin \(g — c)^f + a — 7} 

4- §m sin {(2 — 2m — g)pt- 2^8 4- 7} 

The expression for s is more complex in th’s form than when 
given in terms of the true longitude 0. 
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MoorHs Parallax. 


58. If P be the moon’s mean parallax, and FT the parallax 
at the time 


radius of earth 
distance of ]) 


y = Bu(l —^5*) to the third order, 

-V(iV) 


c= JBm {1 — cos2 {gpt — 7)} to the second order, 

1 1 - — e* + e cos(cpi — a) 4- e* C082 (cpt - a) 

4* m? cos {(2 — ^ 7 ri)pi - 2/3} 

4 ^^5wccos{(2--2??t — c)j9^ — 2)84-a} ; 
but P= the portion which is independent of periodical terms, 
= Ra (1 — — e^) ; 

i l 4- e cos(cp^ — a) 4- e* cos2 (cp^ — a) 

4- cos {(2 — 2m) pt — 2/8} 

4- {(2 •“ 2m — c)pf — 2)8 + a} 

neglecting terms of the third order. We see that, to the 
second order, the variable part of the parallax is independent 
of the inclination. 


59. Here we terminate our approximations to the values 
of w, 5 , and 9. If we wished to carry them to the third 
order, it would be necessary to include some terms of the 
fourth and fifth orders according to Art .(30), and the ap- 
proximate values of P, P, and 8^ given in Art. ( 24 ), would 
no longer be sufficiently accurate, but we should have to 
recur to the exact values, and from them obtain terms of 
an order beyond those already employed.* The process 
followed in the preceding pages is a sufficient clue to what 
would have to be done for a higher approximation. 


See Parallactic InequalUyy Art. (106). 
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The coordinates w' and & of the sun’s position are, by 
the theory of elliptic motion, known in terms of the time % 
and t is given in terms of the longitude d by the equation 
02. Hence w' and & can be obtained in terms of 0} but 
it will be necessary to take into account the slow progressive 
motion of the sun’s perigee, which we have hitherto neglected. 
This we may do by writing c0' — f for 0' — c* being a 
quantity which differs very little from unity.* 

These values of u\ 0\ together with those of u and $ 
given by and 8^^ are then to be substituted in the 
corrected values of the forces, and thence in the differential 
equations. The integrations being performed as before will 
give the values of 5, and t in terms of 0 to the third 
order, and from these, as in Arts. (55), (56), and (57), may bo- 
obtained «, 5, and 0 in terms of U 

60. More approximate values of c and g are obtained at 
the same time, by means of the coeflGicients of cos (c0 — a} 
and sin (g0 — 7) in the differential equations, (see Appendix,, 
Arts. 102 and 103). 

61. The values to the fourth order are then obtained 
from those to the third by continuing the same process, and 


* ‘En reflechissant sur lestermes quedoivent introduire toutes Ics quantitca 

* pr^cedentes, on voit qu’il se pent glisser des cosinua de Tangle 0 dont nous avona 
*yvL le dangcreux effet d’amener dans la valeur de u des arcs au Teu de Icura 

* cosinus j de tcls termes viendront, par exemple, de la combinaison des cosinus 

'de (1 — m) 0 avec des cosinus de md 

* Pour eviter cet inconvenient qui oterait ^ la solution precedento 

'Tavantage de convenir k un aussi grand nombre de revolutions qu*’on voudrait, 
' et la priverait de la simplicitc et de T universalite si precieuses en matbematiques, 
'il faut commencer par en chercher la cause. Or, on decouvre facilement que 
'ces termes ne viennent que de ce qu'on a suppose fixe Tapogee du soleil, ce 

* qui n’est pas permis en toute rigueur, puisque quelque petite que soit sur cet astro 
‘ Taction de la lune, elle n’en est pas moins reelle et doit loi produire un mouvement 

* d’apogee quoique tri)s lent k la verity’.— Clairaut, Theorie de la Lune^ p, 6d. 
2me Edition, 
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SO on to the fifth and higher orders ; but the calculations are 
so complex that the approximations have not been carried 
beyond the fifth order, and already the value of 6 in terms 
of t contains 128 periodical terms, without including those 
due to the disturbances produced by the planets. The coeffi- 
cients of these periodical terms are functions of w, e\ ^ , 

c, g, /c, and are themselves very complicated under their 
literal forms : that of the term whose argument is twice the 
difference of the longitude of the sun and moon, for instance, 
is itself composed of 46 terms, combinations of the preceding 
constants. 

See Pontdcoulant, Syst^me du Monde^ tom. iv., p. 572. 
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CHAPTER V. 

NUMERICAL VALUES OF THE COEFFICIENTS. 

62. Having thus, from theory, obtained the form of the 
developments of the coordinates of the moon’s position at 
any time, the next necessary step is the determination of the 
numerical values of the coefficients of the several terms. 

We here give three different methods which may be 
employed for that purpose, and these may, moreover, be 
combined according to circumstances. 

63. First Method, The values of the constants tw, a, 

yS, 7 , f which enter into the aryuments^ and of the additional 
ones a, 7c, &c. which enter into the coefficients^ of the 

terras in the previous developments, may be obtained with 
great accuracy from observation. 

For this we must employ observations separated by very 
long intervals, such, for instance, as ancient and modern 
eclipses; also particular observations of the sun and moon 
made when the bodies occupy certain selected positions; 
and other observations of a special character. 

When the values of the elements have been so obtained, 
then the theoretical values of the coefficients may be com- 
puted by substitution in the analytical expressions. 

64. Second Method, Let the constants which enter into 
the arguments be determined as in the first method ; and let 
a large number of observations be made, from each of which 
a value of the true longitude, latitude, or parallax is obtaine^^ 
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togctlier with the corresponding value of t reckoned from 
the fixed epoch when the mean longitude is zero. Let these 
corresponding values be substituted In the equations, each 
observation thus giving rise to a relation between the unknown 
constant coefficients, 

A very great number of equations having been thus 
obtained, let them, by the method of least squares or some 
analogous process, be reduced to as many as there are 
coeflScients to be determined. The solution of these simple 
equations will give the required values. 

This method, however, would scarcely be practicable in 
a high order of approximation. For instance, in the fifth 
order, as stated in Art. (61), each observation would give 
rise to an equation containing 129 unknown quantities and 
the immense number of equations so obtained would have 
to be reduced to 129 equations of 130 terms each. 

65. Third Method, When the constants which enter into 
the arguments have been determined by the first method, we 
may obtain any one of the coeflScients independently of all 
the others by the following process, provided the number of 
observations be very great. 

Let the form of the function be 

F= A + ^sin^+ (7sin0 + &c-) 
and let it be required to determine the constants (7, &c. 

separately; <^, Ac., being known functions of the time. 

Let the results of a great number of observations corre- 
spoiiding to values 6^^ 6',, «&c., «&c., be F„ F^y 

&c. ; so that 

= J + J?sln0j+ (7sin<^j-f &c., 

F2 = A + J5sln02+ f7sin^2 + &c., 

F3 = A + £sin 03 + (7sln<^3 + &c., 


F„ = A -f i? sin 6^^ + C sln<f>^ + &c. 
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Now, n being very great^ we may assume that the sum of 
the positive values of each periodical term will be about 
counterbalanced by the sum of its negative values; and 
therefore, that if we add all the equations together these 
terms will disappear ; 


therefore 


VAV,+y,+ +v. 


which determines the non-periodic part of the function. 

To determine B, Let the observations be divided into 
two sets separating the positive and negative values of sin 6 ; 
then the other periodical terms, not having the same period, 
may be considered as cancelling themselves in adding up the 
terms of each set. Let there be r terms in the first set 

and s terras In the second, and let P, P', V'" be the 

values of V corresponding to positive values of sin0, which 
values we may assume to be uniformly distributed from 

sin^ to sinTT, and therefore to be sin 80, sin 280, sin r. 80, 

where r.80 = 7r. 

And, again, let F^, F^,, F^,,, F,, be the values of F 

corresponding to the negative values of sin0 ; viz., sin (- A0), 
sin (- 2 A0), sin (- 5. A0), where 5. A0 = tt. Then, 


F' =-4 + J5sln80 -f C^sIn</>' F^ =^4- J^sln A0 + (7sin(^^ 4..., 
P' = A + 58in2.804 (Jslnf'+..., F,^ = il-i?sin2. A04 C'sin0,,4..., 


F*" =A + i?sInr.80+ Csin^*' V, = A-jBsins.A0+ Csin^, 


therefore 


therefore 


p+ P'+...+P=r.A+J5s7(sIn0) F, + F, ,+...+ F,=5.A-J8S/(sm0) 


= rA+gjV( 8 m 0 , 80 ); 




therefore 


therefore 


r irj„ s irj„ 


. 22 ? 

=-4+ — . 

IT 


-A — ■ 

IT 
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therefore B = 


z (Xl±Il±ii±I: 

4 V r 


s )' 


and in a similar manner may each of the coeflScients be 
independently determined.* 


66. When the periods of two of the terms differ bu^ 
slightly — for instance if 0 and (f) go through their periodic 
variations very nearly in the same time, — the method could 
not then with safety be applied ; for, since the same values 
of 6 and (f) would very nearly recur together during a longer 
time than that through which the observations would ex- 
tend, the two terms would be so blended in the value of 
V that they would enter nearly as one term — the difference 
between 0 and ^ would be very nearly the same at the 
end as at the beginning of the scries of observations. 


If r and 8 are not Bufliciently great to allow us to substitute ein Odd for 
sin 6 . SOf we must proceed as follows : 

V' 4- V" 4- + V*' — rA + B (sin 4- sin 250 4- 4- Binr50) 

^ sin ^ (r 4- 1 ) 50 sin ^7' , 50 
Bin>i50 


= rA-\-B~ 


r'4- V" 4-.. 


, , „cos|50 
-= rA 4- B 'V-, , 

Bin ^50 

Vr B 1 - 1^502 

= ^ + ~ i-Tm > nearly, 


" A + 


= A 4- 




Similarly, 

therefore 


r 4- V, 4- 4- Vs 




B 

TT 

/F'4-F"+... 

,,.4- V, 4- V,, + 4- r 

4 1 


rv 

s 

'M 

24 ((•“■‘‘sV. 
r i\ i\ 

\ 

1 fV’+V"+... 

...4- F'' F, 4-F,, 4- 4- 

“4 1 


8 
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67. Let us suppose the periods to be actually identical, 
so that ^ = 0 4- a, a being some constant angle ; then 

B&ind-hG sin (j> 

may be written {B+C cos a) sin ^ + C sin a cos 

F=-4 + {i?+ (7cosa) sln^+ (7 sin a cos 

If now we divide the observations, as before, into two 
sets, corresponding to the positive and negative values of 
sin0, the terms involving cos0 will disappear in the sum- 
mation of each set ; and, following the process of the method, 
we shall find the value of 


B-\- G cos a = J/ suppose. 

Dividing again into two sets corresponding to the positive 
and negative valuta of cos0, the terms in sin 0 will bo can- 
celled, and the same process will give the value of 

(7 sin a = iV suppose. 

Treating the observations in the same way with respect 
to the angle <j?), wc get two results, 

(7+i?cosa = if', 

- B sina — N’ • 
from these four equations we easily get 


= -TTv TV or = iV , 


Mj Nj N' arc connected by the equation of condition, 


68. When the periods of d and ^ arc nearly, but not 
exactly, the same, this equation of condition will not hold, 
and the preceding values of B and C would not be exactly 
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correct ; but yet they would be very approximate, especially 
if the mean between the two values of B be taken. 

Or, we may, after having taken one of these slightly 
erroneous values for J?, make a further correction by estab- 
lishing as it were a counterbalancing error in the value of 
(7. Let jy be the value so found for .B; then, from the 
V of each of the observations subtr^-ct the value B' sin 0^ 
the result U will be very nearly equal to ^4 + (7 sin ^ + &c., 
and from the n equations 

— A + (7sln<^j+ 

= (7sin0jj-l- 

(7sin^^+ 

a value (7 of (7 will be obtained, by the rule, which will bo 
very approximate, and, at the same time, agree better with 
B' in satisfying the equations than G itself would do. 

69. When two terms whose periods are nearly equal do 
occur, it is plain, by examining the values of M and M\ that 
the errors which would be committed by following the rule, 
without taking account of this peculiarity, would be the taking 
B + <7 cos a and (7+ B cos a for B and C respectively. 
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CHAPTER VI. 

THYSICAL INTERPRETATION. 

70. The solution of the problem which is the object of 
the Lunar Theory may now be considered as effected ; that 
is, wc have obtained equations which enable us to assign 
the moon’s position in the heavens at any given timo to 
the second order of approximation ; we have explained how 
the numerical values of the coefficients in these equations 
may be determined from observation; and we have, more- 
over, shewn how to proceed in order to obtain a higher 
approximation.* 

It will, however, be interesting to discuss the results wc 
have arrived at : — to see whether they will enable us to form 
some idea of the nature of the moon’s complex motion ; also 
how far they will explain those inequalities or departures 
from uniform circular motion which ancient astronomers had 
observed; but which, until the time of Newton, were so 
many unconnected phenomena; or, at least, had only such 
arbitrary connexions as the astronomers chose to asrign, by 
grafting one eccentric or epicycle on another as each newly 
discovered inequality seemed to render it necessary. 

It is true that our expressions, composed of periodic 
terms, are nothing more than translations into analytical 


* The means of taking into account the ellipsoidal figure of the earth and the 
disturbances produced by the planets, are too complex t > form part of an intro- 
ductory treatise. For information on these points reference may be made to 
Airy'a Figure of the Earth, PontecoulanPs Sy$tme <iu Monde^ vol. iv. 
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language of the epicycles of the ancient;* but they arc 
evolved directly from the fundamental laws of force and 
motion, and as many new terras as wc please may be ob- 
tained by carrying on the same process ; whereas the epicycles 
of Hipparchus and his followers were the result of numerous 
and laborious observations and comparisons of observations , 
each epicycle being introduced to correct Its predecessor 
when this one was found inadequate to give the position 
of the body at all times: just as with us, the terms of the 
second order correct the rough results given by those of 
the first; the terms of the third order correct those of the 
second, and so on. Hut it is impossible to conceiv^e that 
observation alone could have detected (dl those minute 
irregularities which theory makes known to us in the terms 
of the third and higher orders, even supposing our instruments 
far more perfect than they are ; and it will filways be a 
subject of admiration and surprise, that Tycho, Kepler, and 
their predecessors should have been able to fvcl their way 
so far among the Lunar inequalities, with the means of 
observation they possessed, 

LONGITUDE OF THE MOON. 

71. We shall first discuss the expression for the moon's 
longitude, as found Art. (55), 

0 2e sin (cpt - a) + sin 2 {ept - a) 

+ Y sin {(2 — 2in - c) pt - 2/3 -f a\ 

+ {(2 - — 2 / 8 j 

- sin {mpt + ~ f) 

- H’’ xin 2 {ffpt- y). 

The mean value of 0 is i^ud in order to judge of the 
elFect of any of the small terms, we may consider them 


* IScc WhcwcU's l/hlnry of thi Indvcfifi: J^^deuas. 
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one »Tt a time as a correction on this mean value or 
we may select a combination of two or more to form this 
correction. 

We shall have instances of combination in the elliptic 
inequality and the eveUion^ Ari^. (73) and (77); but in the 
''eniaining inequalities each term of the expression will form 
a correction to be considered by itself. 

72. Neglecting all the periodical terms, we have 

(10 

dt 

which indicates uniform angular velocity; and as, to the 
same order, the value of n is constant, the two togetlier 
indicate that the moon moves uniformly in a circle, the 

27r 

period of a revolution being - , which is, therefore, the 

expression for a mean sidereal month, or about 27^^ days.^* 
The value of^; is, according to Art. (54), given by 

and as in is due to the disturbing action of the sun, we sec 
that the mean angular velocity is less, and therefore the mean 
periodic time greater than if there were no disturbance. 

Elliptic Inequality or Equation of the Centre. 

73. We shall next consider the eflfect of the first three 
terras together : the effect of the second alone, as a correction 

pt^ will be discussed in the Historical Chapter, Art. (115). 


♦ The accurate value was 27d. 7h. 43m. ll*2GIs. in the year 1801. See 
Art. (108). 


K 
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6 -j)t + 2c sin (cpt ~ a) + \e' sin 2 (cpt - a), 
which may be written 

6=zpt-\-2e sin {a -f-(l — sin 2{pt — {«+ (1 —c)pt]'\. 

But the connexion between the longtltudc and the time in 
an ellipse described about a centre of force in the focus i’’^ 
Art. (13), to the second order of small quantities: 

0 = n# -f 2e sin (nt - a') + sin 2 (nt — a'), 

where n is the mean motion, e the eccentricity, and a the 
longitude of the apse.^ 

Hence, the terms we arc now considering indicate motion 
in fin ellipse ; the mean motion being p, the eccentricity c, 
and the longitude of the apse a4(l — c) that is, the 
apse is not stationary but has a progressl re motion in longi- 
tude, uniform, and equal to (1 ~ c) p, 

74, The two terms 2c sin (cpt — a) + f sin 2 (cpt — a) con- 
stitute the elliptic inequality^ and their effect may be further 
illustrated by means of a diagram. 

Let the full line A MB re- 
present the moon’s orbit about 
the earth when the time t 
commences, that is, when the 
moon’s mean place is in the prime 
radius from which the lon- 
gitudes are reckoned. 

The angle *TEA^ the longi- 
tude of the apse, is then a. At the time when the moon’s 
mean longitude Is TEM^pt^ the apse line will have moved 

♦ The epoch « which appears in the expression of Art. (13) is here omitted ; 
a proper assumption for tho origin of f, as explained in Art. (36), enabling us 
to avoid the g. 
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in the same direction through the angle AJEA' = (1 — c) TJSilfj 
and the orbit will have taken the position indicated by the 
dotted ellipse. The true place of the moon in this orbit, 
so far as these two terms are concerned, will be ?«, where 

MUm = 2e sin (cpt — a) + f 6* sin 2 (cpt — a) 

= 2(3 sit A' EM 4 mi2AEM 
= 2e sin-4'^Jf (1 + cos AE3I } ; 

which, since e is about is positive from perigee to apogee, 
and therefore the true place before the mean; and the 
contrary from apogee to perigee: at the apses the places 
will coincide. 

75. The angular velocity of the apse is (1 — c) p, or, if 
for c we put the value found in Art. (52), the velocity will 
be Hence, while the moon describes 360'', the apse 

should describe |m^3G0'’= If’ nearly, rn being about 

But Hipparchus had found, and all modern observations 
confirm his result, that the motion of the apse is about S'' in 
each revolution of the moon. See Art. (IIB). 

This difference arises from our value of c not being repre- 
sented with suflScient accuracy by 1 — \rn\ 

Newton himself was aware of this apparent discrepancy 
between his theory and observation ; and we are led, by his 
own expressions (Scholium to Prop. 85, lib. iii. in the first 
edition of the Principia)^ to conclude that he had got over 
the diflSculty. This is rendered highly probable when we 
consider that he had solved a somewhat similar problem 
in the case of the node ; but he has nowhere given a state- 
ment of his method : and Clairaut, to whom we are indebted 
for the solution, was on the point of publishing a new hy- 
pothesis of the laws of attraction, in order to account for 
it, when it occurred to him to carry the approximations ^o 
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the third order, and he there found the next term in the 
value of c nearly as considerable as the one already obtained. 
See Appendix. The value is 

c = 1 ~ • 

I - c = = f wi® (1 + 5 

(1 ~ c) 3G{/ = (1 4- iW) {value found previously) 

= 22" nearly, 

thus reconciling theory and observation, and removing what 
had proved a great stumbling-block in the way of all as- 
tronomers.* 

When the value of c is carried to higher orders of ap- 
proximation, the most perfect agreement is obtained. 

The motion of the apse line is coiibldercd by Newton 
In his rrinclj^ia^ lib. I., Prop. 6G, Cor. 7. 


Evection. 

76. The next term + sin {(2 — 2m - o) pt - 2/3 + a} 
in the value of 6 has been named the Evection, We shall 
consider its effect in two different ways. 

First, by itself, as forming a correction on 

6 —pt 4- sin {(2 - 2m — c) pi — 2/3 4- «}• 

Let ]) = pi = moon’s mean longitude at time 

^ = mpt 4- i8 = sun’s 

a = (1 - e) a = mean longitude of apse , 

then 

d~pt-V sin [2 { pt - {mpt 4- /3)} — [pt - (1 - c) pi 4- a}] 
-pt 4- sin [2 (j - 0) ~ (5 - a)}. 


>^00 Dr. Wlicwdl’s? Brifl^cwaUr Tnatin', 
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The eflfcct of this term will therefore be as follows : 

In sjzygies 

0 :=ipt — \^m€ sin (]) - a ) ; 

or the true place of the moon will be before or behind the 
mean, according as the moon, at the same time, is between 
apogee and perigee or between perigee and apogee. 

In quadratures 

0 -f sin (J — a'), 

and the circumstances will be exactly reversed. 

In both cases, the correction will vanish when the apse 
happens to be in syzygy or quadrature with regard to the 
sun at the same time as the moon. 

In intermediate positions, the nature of the correction 
is more complex, but it will always vanish when the sun 
is at the middle point between the moon and the apse, 
or when distant 90^ or 180"* from that point; for if 

G = “ r . 90°, where r = 0, ± 1, or 2, 

sin [2 (5 - g) ~ (3) - a)] = sin (]) + a' - 2o) 

= sin r. 180° 

= 0 . 

77. The other and more usual method of considering 
the effect of this terra is in combination with the two terms 
of the elliptic inequality, as follows : 

To determine the change in the position of the apse and in 
the eccentricity of the moon^s orbit produced by the emotion. 

Taking the elliptic inequality and the cvection together, 
we have 

0^pt’\- 2e sin - a) + ^e^ sin 2 {cpt - a) 

+ ^me sin {{2 — 2m - c)pt — 2^-\‘ a]. 
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Let a' be the longitude of the apse at time t on supposition 
of uniform progression, 

O sun ; 

whence a = (1 —c)iA + a, 

O = milt + yS. 

And the above may be written ^ 

6 —pt -f 2e sin [cpt — a) -f sin 2 {cpt — a) 

+ sin [cpt - a 4- 2 (a - g)} ; 
and the second and fourth terms may be combined into one, 
2E sin [cpt — a 4 S), 

if we assume cos 3 = e 4 cos 2 (a ~ G ), 

and E sin 8 = sin 2 (a' - g) ; 

whence tan 8 and E may be found ; and approximately, 

E=e {1 4 V'wi C032 (a' — g)1) 

8 = sin 2 (a — g). 

The term \e? fkii2[cpt-‘oi) may also, therefore, to the 
second order, be expressed by 

^E'^ sin 2 [cpt — a 4 8) , 
and the longitude becomes 

^ 4 2E sin [cpt ~ a 4 8) 4 %E'^ sin 2 [cpt - a 4 8), 

or d^pt 4 2i? sin [iit - a 4 8) 4 \E'^ sin 2 - a 4 8) ; 

but the last two terms constitute elliptic inequality in an orbit 
whose eccentricity is E and longitude of the apse a' — 8; 
therefore the evcction, taken in conjunction with elliptic 
inequality, has the effect of rendering the eccentricity of the 
moon’s orbit variable, increasing it by when the apse- 
line is in syzygy, and diminishing it by the same quantity 
when the apse-line is in quadrature ; the general expression 
for the increment being cos 2 (a — g). 
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And another effect of this terra is, to dirainisli the longitude 
of the apse, calculated on the supposition of its uniform 
progression, by the quantity 8 = sin2 (a - O) ; so that 
the apse is behind its mean place when in the first and third 
quadrant in advance of the sun, and before its mean place 
in the second and fourth. 

The cycle of these jhanges will evidently be completed 
in the period of half a revolution (if the sun with respect 
to the apse, or in about -/j. of a ycar.*^’ 

78. The period of the evcction itself, considered inde- 
pendently of its effect on the orbit, is tlie time in which 
the argument (2 — 2m — c)])t — 2y8 + a will increase by 27r. 

Therefore period of evection 

27r mean sidereal month 

(2 — 2m — c) p 2 — 2in — c 

mean sidereal month _ 27 days 
i — 2iii + 1 ■“ A ’ 

= 31^ days, nearly, f 

Newton has considered the evection, so far as it ari^.cs 
from the central disturbing force, in Prop. G6, Cor, 9, of 
the Fnncijna. 


nearly, 


Variation, 

79. To explain the physical rneaniwj of the term 
y w* sin {(2 — 2m) pi •— 2yQ}, 

* The change of eccentricity and Uie variation in tlie motion of the apso 
follow the same law as the abscif-sa and ordinate of an ellipse referred to 
its centre ; for if E — x and o-y^ then 

t Tlie accurate value is ill ’8119 days. 
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in the expression for the moon'^s longitude^ 

0 ^pt + y sin {(2 — 2m) pt — 2^}. 

Let }) represent the moon’s mean longitude at time tj 

0 sun's 

therefore }) —pt^ 

O = mpt + 0 ; 

and the value of 0 becomes 

q), 

which shews that from syzygy to quadrature the moon’s 
true place is before the mean, and behind it from quad- 
rature to syzygy; the maximum difference being in 

the octants. 

The angular velocity of the moon, so far as this term 
is concerned, is 

=:p -f Y (1 - m) vi^p cos2 (]) - O), 

=/? {1 + cos2 (]) - g)}, nearly, 

which exceeds the mean angular velocity p at syzygics, is 
equal to it in the octants, and less in the quadratures. 

This inequality has been called the Variation^ its period 
is the time in which the argument (2 — — 2^ will 
increase by 27r ; 

. , ^ - 27r mean synodical nmnth 

••• — 2 

= 14| days, nearly.* 

80. The quantity y wi’ is only the first term of an endless 
scries which constitutes the cocflicient of the variation, the 
other terms being obtained by carrying the approximation 


The accurate value is 11*765294 days. 
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to a higher order. It is then found that the next term in 
the coefficient is f ) which is about y\- of the first term ; 
and as there are several other important terms, it is only 
by carrying the approximation to a higher order (the 5th 
at least) that the value of this coefficient can be obtained 
with sufficient accuracy from theory. In fact, V«i’ would 
give a coefficient of 20' 27" only; whereas the accurate value 
is found to be 

The same remark applies also to the coefficients of all 
the other terms. 


81. As far as terms of the second order, the coefficient 
of the variation is independent of e the eccentricity, and k 
the inclination of the orbit. It would therefore be the same 
in an orbit originally circular, whose plane coincided with 
the plane of the ecliptic: it is thus that Newton has con- 
sidered it. Princlp. prop. GO, cor. 3, 4, and 5, 


Annual Jupiation, 


82. To explain the pliysiml meaning of the term 
— 3m// sin [nipt -f /3 - f) 
in the exjiression for the 'nwoii's longitude, 

— *6me sin [mpt -f yS - f ), 

—pt — 3m/^' sin (longitude of sun — longitude of sun’s perigee), 
^pt - onie sin (sun’s anomaly). 

Hence, while the sun moves from his perigee to his 
apogee, the true place of the moon will be behind the mean ; 
and from apogee to perigee, before it. The period being an 
anomalistic year, the effect is called Annual Eguation, 
Differentiating 6 we get 


de 

dt 


—p {1 — 3//iV cos (sun’s anomaly)}. 
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Hence, so far as this Inequality Is concerned, the moon’s 
angular velocity Is least when the sun Is in perigee,' that 
is at present about the 1st of January, and greatest when 
the sun is in apogee, or about the Ist of July. 

The annual equation is, to this order, Independent of the 
eccentricity and inclination of the moon’s orbit, and therefore, 
like the variation, would be the same in an orbit originally 
circular. Vide Newton, IVincipia, prop. 66, cor. 6. 

Reduction. 

83. Before considering the effect of the term 

A** 

- jBm2{(jpt-ry)^ 

which, as wo shall see Art. (89), is very nearly equal to 
the difference between the longitude in the orbit and the 
longitude in the ecliptic, it will be convenient to examine 
the expression for the latitude of the moon, and to see how 
the motion of the node is connected with the value of g. 

LATITUDE OF THE MOON. 

84. The expression found for the tangent of the latitude,* 
Art. (51), is 

s = A: sin [gO - 7) 4 Imk sin {(2 - 2Tn - <7) ^ - 2/8 + 7}. 

If wo reject all small terms, we have 

a = 0, 

or the orbit of the moon coinciding with the ecliptic, which 
is a first rough approximation to its true position. 


♦ This expression for the tangent of the latitnde is more conrenient than 
that which gives it in terms of the mean longitude, Art. (57), on account of 
the 'css number of terms involved. See Pontecoulart, vol. iv. p. 030. 
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85. Taking the first term of the expansion 
5 = A sin (^^ — 7 ), 

we may write it 

« = A;8:n 

Let TNm be the ecliptic, N the 
moon’s node when her true longitude 
is aero, and let M be the position 
of the moon at time ty in her place 
referred to the ecliptic ; 

therefore T-N'** 7 , = tanil///e = j?. 

Take NN' — {g-l)0 in a retrograde direction; and join 
MN' by an arc of great circle ; 

then tan -Min = tani/-^'msinj^'?}i, 
or s = tan MN'm sin [0 — {7 — - 1 ) 0 ]] ; 

which, compared with the value of s given above, shews that 
MN'm = tan'^A is constai.t, and therefore the term k sin {fj0—y) 
indicates that the moon moves in an orbit inclined at an 
angle tan”^^ to the ecliptic, and whose node regredes along 

the ecliptic with the velocity with a mfian 

velocity {g - 1 )^. 



86. Hence the period of a revolution of the nodes 
27r one sidereal month 

«/-i ’ 

but, from Art. (51), the value of = 1 + ; 


_ one sidereal month 

= 6511 days, nearly. 

This will, for the same reason as in the case of the apse. 
Art. (75), be modified when we carry the approximation 


therefore period of revolution of nodes 
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to a higher degree ; tliis value of cj Is, however, much more 
accurate than the corresponding value of c, for the third 
term of cj is small ; the value to the third order being (see 
Appendix, Art. 102) 

, , . 1 It I . I. , 1 one sidereal month 

and the period or revolution ot the noc^es o — 

= G7()5 days, nearly. 

This is not far from the accurate value as given by 
observation, and when the approximation to the value of g 
is carried to a higher order, the {igreemcnt is nearly perfect. 

The true value is 6793*39 days, that is about ISyrs. 7 mo. 

Ecectioa in Latitude. 

87. To explain the variation of the inclination and the 
irregularifg in the ^notion of the node expressed hj the term 
+ ^nih sin {(2 - 2ni — y) 0 — 2/S + 7]. 

This term, as a correction on the preceding, is analogous 
to the evection as a correction on the elliptic inequality. 

Taking the two terms together, 

s = k sin {gd — 7) + %ink siii {(2 — 2m — 0 — 2)8 + 7}. 

Let ]) = longitude of moon = 6^ 


0= sun =7R^ + )8, 

node =7 — (<7-l)0; 


therefore s — k sin (5 ~ + Inik sin [I) — — 2(0 — SI))}* 

Now these two terms may be combined into one, 

,s* = K sin (]) - Q - A), 
if K cos A = k + ^mk cos2 (O — g^), 

K sin A == ^mk sin 2 (O - Q), 
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whence A and K may be found ; and, approximately, 

4- co92(0- S^)}, 

A = gw^sill2(0 - Q), 

but the equation 

.9 = JTsinC})- Q - A) 

represents motion in ai» orbit inclined at an angle tan'^JT to. 
the ecliptic, and the longitude of whose node is Q + A. 

This term has therefore the following effects : 

1st. The inclination of the moon’s orbit is variable, its 
tangent increases by when the nodes are in syzygies, 
relatively to the sun, and decreases by the same quantity 
when they are in quadrature ; the general expression for 
the increase being ‘^mk cos 2 (^ 1 ’ — Q). 

2nd. Tlie longiUide of the node, calculated on supposition 
of a uniform regression, is increased by A = sin2 
so that the node is before its mctan place while in the first 
or third quadrant In front of the sun, and behind it in the 
second and fourth. Pruicipla^ book III., props. 33 and 35. 

The cycle of these changes will be completed in the 
period of half a revolution of the sun with respect to the 
node, that is, in 173*21 days, not quite half-a-year. 

88. The tangent of the latitude has here been obtained ; 
if we wish to have the latitude itself it will bo given by the 
formula 

latitude = s — J.s’* + -}/' - &c, 

which, to the degree of approximation adopted, will clearly 
be the same as a. 


lieduction, 

80. We may now consider the term which we had 
neglected (Art. 83J in the expression for the longitude, namely, 
— sin 2 {(jpt — 7). 
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Let N be the position of the node when the moon^s 
longitude is <?, M the place of the 
moon, m the place referred to the 
ecliptic. 

Therefore Tm = 

N‘m = ““ 7 , 

tan^= A. 



The right-angled spherical triangle N3{jn gives 

trin Nin ^ 

VAtiNJir 

tan NM— tan Nin 


therefore 


or 


cos N= 
1 - cosN 


1 + cos N tan NM + tan Nm ’ 

j, -AT ^ sin (NM — Niii) 

2 sin (NM + Nm) ’ 


or, since both N and NM—Nni are small, 


tan^iY NM- Nm . , , 

“4 “ 


therefore NM-Nm=lk^ sm2(g0—y)=lk^ sin2({jpt—y)y nearly. 

Hence this term, which is called the reduction^ is approxi- 
mately the dilFerence between the longitude in the orbit and 
the longitude in the ecliptic. 


RiVDIUS VECTOR. 

90. To explain the physical meaning of the terms in the 
value of u. 

Wo shall, for the explanation, make use of the formula 
which gives the value of u in terms of the true longitude, 
Art. (52). 
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Neglecting the periodical terms, we have for the mean 
value 

= a (I — f 

The term which is a consequence of the disturbing 

effect of the sun, shews that the mean value of the moon’s 
radius vector, and tliereforc the orbit itself, is larger than if 
there were no disturbance. 

Elliptic Inequality, 

91. To explain the effect of the term of the first order^ 

u^a[l + 6CO8(c0 — a)}, 

= a [1 cos^— {a4- (l - c) 6^jj. 

This is the elliptic inequality^ and indicates motion in an 
ellipse whose eccentricity is e and longitude of the apso 
a + (l— c)0; and the same conclusion is drawn with respect 
to the motion of the apse as in Art. (73). 

Evcction, 

92. To explain the physical meaning of the term 

^mea cos {(2 — 2m — c) 0 — 2^ + a}. 

This, as in the case of the corresponding term in the longi- 
tude, is best considered in connection with the elliptic 
inequality, and exactly the same results will follow. 

Thus, calling 5, G), and a' the true longitudes of the moon, 
sun, and apse (the last calculated on supposition of uniform 
motion), these two terms may be written, 

w = a[l + eco3(}) -a')— ^me cos {D — a' -f 2 (a' — O))] 

= a [1 + j& cos (D — a' + 8)] ; 
where E cos 8 = e -f V 2 (a — g)> 

jS sin 8 = i/me sin 2 (a - O). 

These are Identical *with the equations of Art. (77). 



so 


LUNAR THEORY. 


Variatioiu 

93. To explain tlia effect of the fenti in\i cos \{2—2ni) ^-2)8}, 

t4 = [I + ni^ cos }(2 - 2m )6 — 2/3}] 

= a [1 -H cos2 (]) - q)J. 

As far as this term is concerned,* the moon’s orbit would 
be an oval having its longest dijimctcr in quadratures and 
least in syzygies. Prindpia^ lib. 1. prop. GG, cor. 4. 

The ratio of the axes of the oval orbit will be 

1 “f" 7n^ 

I ‘2= wll nearly, m being ‘0748. 

See Prineijna^ lib. lir. prop. 28. 

Jledurtion. 

94. The last important periodical term in the value of u is 

/ /I ^ 

— j cos 2 (</£/ — 7). 

This term expresses approximately the variation in the 
diiltTence between the values of u in the orbit and in the 
ecliptic. 

For if ?q be the reciprocal of the value of the radius 
vector in the orbit, 

7q = u cos (latitude), 

= V = (1 — nearly ; 

VU + .s'-j ^ 

therefore u — w, = J iis^ == siu^ (r/6 — 7) 

= cos 2 {gd - 7) 

= const. — cos 2 {gd — 7). 

Taking this result in connexion wdth that of Art. (89) 
we sec that the reduction^ so far as periodical terms to the 
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second order arc concerned, is simply a geometrical con- 
sequence of the inclination of the orbit; and that, if we 
measured the longitude and the radius vector along the 
orbit instead of taking their projection on the piano of the 
ecliptic, these periodical terms would not appear. 

95. There are no ''thcr terms of the second order In the 
value of The annual equation^ which, in the longitude, 
is of the second order, is only of the third order in the 
radius vector. 


Periodic time of the Moon, 

9G. We have seen. Art. (72), that the periodic time of 
the moon is greater than if there w(‘,rc no disturbing force ; 
but this refers to the mean periodic time established on an 
interval of a great number of years, so that the circular 
functions in the expression are then extremely small com- 
pared with the quantity pt which has uniformly increased. 

When, however, we consider only a few revolutions, these 
terms may not all be neglected. The eMiptio inequaliUj and 
the evcction go through their values in about a month,* the 
variation and reduction in about half-a-month ; their effects, 
therefore, on the length of the period can scarcely be con- 
sidered, as they will increase one portion and then decrease 
another of the same mouth. 

But the annual equation takes one year to go through its 
cycle, and, during this time, the moon has described thirteen 
revolutions; hence, fluctuations may, and, as we shall now 
shew, do take place in the lengths of the sidereal months 
during the year. 

We have, considering only the annual equation. Art. (82), 
pt ?,me mi(md + 0 ' - 

M 
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Let T be the length of the period ; then, when 6 is 
increased by 27r, t becomes T] 
therefore (<+?') = 27r + 0 + Zme sin (2m7r + mO + /3 - 
whence 27r + Gme sin imr cos (inir + 7116 + — 0 5 

therefore T=^ mean period + sinwiTr cos(o — Oj 

^ ( 

where o + + W27r = sun’s longitude at the beginning 

of the month + ??^7^ 

= sun’s longitude at the middle of the month. 

Hence T will be longest when 0 — f=0, 
and shortest when O — f = tt ; 

or T will be longest when the sun at the middle of the month 
is in perig(jc, and shortest when in apogee ; but, at presei^t, 
the sun is in perigee about the 1st of January, and in apogee 
about the 1st of July; tlicrcforc, owing to annual equation^ 
the w'inter months will be longer than the summer months, 
the difference between a sidereal moiuh in January and July, 
from this cause, being about 20 minutes. 

97. All the inequalities or equations, which onr expres- 
sions contain, have tlius received a physical interpretation. 
They were the only ones known before Newton had estab- 
lished his theory, but the necessity for such corrections was 
fully recognized, and the values of the coefficients had already 
been pretty accurately determined ; still, with the exception 
of tlic reduction, which is geometrically necessary, they were 
corrections empirically made, and it was scarcely to be ex- 
}'cctcd that any but the larger inequalities, viz. those of the 
first and second orders which we have here discussed, could 
be dete(‘tcd by observation : w^e find, however, that three 
others, have, since Ncwtoirs time, been indicated by obser- 
vation before theory had explained their cause. These arc — 
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the secular acceleration^ discovered by Halley ; an inequality, 
found by Mayer, in the longitude of the moon, and of which 
the longitude of the ascending node is the argiuneiit; and 
finally an inequality discovered by Biirg, which has only of 
late years obtained a solution. For a further account of 
these, as also of some other inequalities which theory has 
made known, sec Appendix, Arts. (108) to (112). 
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APPENDIX. 

In this chapter will be found collected a few propositions 
intimately connected with the results or the processes of the 
Lunar Theory as explained in the previous pages. Refer- 
ence has been made to some of them in the course of the 
work, and the interest and importance of tlie others are 
sufficient to justify their introduction here. 


1)8. The moon is roUiined in her orhit Inj the force of 
gravif}/^ that isj hi/ the same force ichich acts on bodies at 
the surface of the earth. 

The proof of this is merely a numerical verification ; the 
data required from observation are, 

the space fallen through from rest in 1" by bodies at the 


earth’s surface = 16*1 feet, 

the radius of the earth = 4000 miles, 

the pcricJic time of the moon ==27j days. 


the distance of the moon from the earth’s ccntre=G0x4000 miles. 

The force of the earth’s attraction cc — ,r . Therefore, the 

(dist.) ' 

space fallen through in 1'" at distance of moon by a body 

1 6*1 

moving from rest under the earth’s action = feet 

= -00447 feet. 



moon’s orbit concave. 


80 


But tlie moon In one second describes an anf^Ic ~ =a>, 

^ 27J.24.G0" ’ 

during which tlie approach to the eartli 

= CO X 4000 X 5280 (vers, cw) feet 

CO X 4000 X 5280 x 27 r' . 

- - 

= *00418 feet. 

Therefore, the space through which the moon is deflected 
in one second from her straight path, is just the quantity 
through which she would fall towards the earth, supposing 
lier to be subject to the earth’s attraction ; and we may, 
therefore, conclude that she is i*etaincd in her orbit by the 
force of gravity. 

When first Newton, in IGCG, attempted to verify this 
result, he found a difference between the two values equal 
to one-sixth of the less: the reason of Iiis failure was the 
incorrect measures of the earth, which he made use of in 
his computation; and it was not till about 16 years later 
that he was led to the true result, by using the more 
correct value of the earth’s radius obtained by Picart. * 

Principia^ lib. III., prop. 4. 

99. The moon\^ orbit is everywhere concave to the sun. 

Let 8^ Ej and M be the centres of the sun, earth, and 
moon. We ought first to apply 
to each body forces equal and 
opposite to those which act on 
the sun in order to bring him to rest. These forces arc, 
however, so small that we may neglect them, and we shall 
consider the moon to be moving about the sun fixed, and to 
be disturbed by the earth alone. 
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The forces on M arc, therefore, 


m 

sfr 


and 




in MS^ 


.in ME, 


This last must be resolved into two, one in MS^ the other 
perpendicular to it. 

Therefore, the whole central force on the moon in MS 


7/1 E , , 

"" 'SAP 


and the proposition will be proved if we shew that this force 
is always positive. Now, Art. (25), 

.... ^ 2ir,Sa'^ 2ir,SM-^ 

periodic time or sun = - ^ — nearly. 


and 

therefore 

therefore 


therefore 


moon 

7/1 , E-^M 


2*17. EM ^ 


SiP EM^ 

m ^ E 

' EAP'> 


nearly ; 


i/i 

SM 

7)1 


, S3r E E 

■i ^ 1 liy iSy- /<M/“ > 


■EAP 


E 




but this is the value of the central force corresponding to 
co3if=-l, and is therefore its least value, lienee the 
central force always tends to the sun, and the path is always 
concave. 

At new moon the force with which the moon tends to the 
sun is, therefore, greater than that with which she tends to 
the earth: the earth being itself in motion in the same 



CENTRAL DISTITRIUXG FORCE. 


87 


direction, and, at that, instant, with p^reater velocity, will 
easily explain how, notwithstanding this, the moon still 
revolves about it.‘^ 


Central and Tangential Disturhing Forces. 


100. We liavo hitherto considered the effects of the 
central and tangential disturbing forces in combination ; but 
it will be interesting to determine to which of them the 
several inequalities principally owe their existence. 

(1) To determine the elfcet of the centnd disturhlng force* 
JIakc T~{)] 


therefore = 0, 

ilo“ k\c 

P 

or substituting for , from Art. 

h tc ' 


dht 


-f u — a 


dd‘^ 
therefore. 


r 1 — :]/t^ ~ \nv + cosc0^ ~ cos (2 — 2m) 6^ 
= 4 ;l?n^e cos (2 — 2 /m — c) 6^^ 


c « .} ''*■ ' V " — 

I -h cos2g0^ ; 

ri — ^“m*4ccos(c^— C()4iW^cos[(2 — 2/M)0-2i8} 

ti — a\ ;4 fe-nie cos {(2 — 2m — c) 0 — 2jS 4 aj 

I cos2 0/0-7). 


If we compare this with the value of u found Art. (52), vre 
see that the elliptic inequality and the reduction are due 
to tlie central or radial force, as also one-half of the varia- 
tion and three-tenths of the evection. 

It would perhaps be proper to separate the absolute cen- 
tral force from the central disturbing force ; the terms due 
to the latter arc those which contain m ; therefore, the elliptic 


* iluj Aul.hor’s A:>trovoniy. p. 233, 
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inequality and tlio reduction are the effects of the former, 
except that in the elliptic inequality the introduction of c, or 
the motion of the apse, is due to the 'disturbing force. * 

(2) To determine the ejj^ect of the ianfjential disturbing force* 
Let the central disturbing force be zero ; 

Pa 

then a = yii (1 neglecting the inclination, and, 


substituting for and for 2 dd their 

values from Art. (r)2J, the equation in u becomes 


d'^u 

10'^ -b w = } t i] cos (2 — 2ih) 9^ -}- ^ ni^e cos (2 — 2m — c) 6 ^ ; 

whence M = a I ‘ ^ 

We have here the remaining half of the variation and the 
remaining seven-tenths of the evection as the effects of the 
tangential disturbance. Also c=l, or, to, the second order, 
the tangential force has no effect on the motion of the apse. 


101. The separate effects of the central and tangential 
forces in producing the inequalities in the longitude may be 
traced in a similar manner. We shall leave this as an 
exercise for the student. The annual 'equation v/ill be found 
to be due to the central force. The other inequalities will 
be divided as In the radius vector; except as regards the 
variation, of which four-elevenths are due to the central 
and the remaining seven-elevenths to the tangential force. 


To calculate the value of g to the third order. 

102. We must hero make use of the results which the 
approximations to the second order have furnished; but as 
the value of g is determined by that terra of the differential 
equation whose frgument is g9—^^ we need only consider 
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those terms which by their combinations will lead to it with- 
out rising to a higher order than the fourth. 

From Arts. (51), (31), wo obtain 

s=-h sin gQ^ + » mk sin (2 — 2)« — g) 

« (1 + cos(2 - 2w) £1,}, 

sin (2 - 1m) 6^. 

From these wc obtain 

T ds 

Jd^'dd -Im) 0, {cos^0, + Jw co%{;i-1m-g) 0,} 

=..,-3*'>iVi sin/;^„ 
d^s 
d& 


+ s = 0, to the second order ; 


\ /*!/'* 

cW fourth order. 

Therefore the equation In 5 becomes, so far as these terms 
are concerned, 

^/^s• 


dff^ 


^S = {- '5?/lVv + d- ? 

assume s = k s\ng 6 ^^ 

therefore (1 -^f/) = “* 

<7 = 1 + '\vi? — . 


To find the value of c to the third order. 

103. Proceeding in a similar manner from the equations, 
Arts. (52), (31), wc get 

= a { 1 + e cos c0, + cos (2 — 2 m) +• Y' cos (2 — 2m — c) , 

N 
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= a - im* -3 {1 + 3 C 08 (2 - 2m) $,}, 
f »«’ ^4 sin (2 - 2m) 

we obtain 

p 

= a - {1 - 3e cosc0 ,...— cos (2 - 27w — c) 

X (1 + 3 cos (2 -27n) 0,} 

= a 4- {3c + Vy^^wc} cosc0,, 

T 

p—g =— {l-4c cosc0j~4m* cos(2-2m) 6-^^me cos(2-2m-c)0J 

/it 

X sin (2 — 2m) 0^ 

= — §m^ {sin (2 - 2m) 0^ — 2e sin (2 — 2m - c) 0^ — sin 
^^=-Gf{c sinc0j + 2?n* sin(2— 2m)0j-f ^-^me sin (2 — 27;i~c)0j] ; 

therefore =. . .|f m*ae cos c0,^ 

the other term Is of the fifth order, 
f T 

Ij3^d0—-i‘^m^ C08(2~2»/i) ^ j~37w*ccos(2-2m-c)0j~4^^m^c cosc^j, 

id0^ 4 w) =« {^ •••- 32 u® cos (2— 2m) 0j-h\^7n^eco3 (2-2m-c)0^; 

therefore 2 4 = — ^m^ae cosc^j, 

and the equation in becomes 
dSi 

dW + ~ ^ *“ cosc^j 4...} 

= a {1 4 (|m*c4 %'^m*c) cosefi^j 4...}. 

Assume m = a (1 4e cosc^j 4...) J 

therefore uc (1 -■ (?) = {|m*c 4 V/m®c) a, 
whence c -U - J/w* — 
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104, Hence, to the third order of approximation, 

mean motion of apse _ 1 — ^ 4 

mean motion of node — 1 ”” wi® 8 — 3we ^ 

and since m = ^ nearly, we see that the moon’s apse 
progredes nearly twice as fast as the node regredcs. 

In the case of one of Jupiter’s satellites, the periodic 
time round Jupiter is only a few of our days, and the 
periodic time of Jupiter round the sun is 12 of our years, 
therefore m, the ratio of these periods, is very small. 

Hence, the apse of one of Jupiter’s satellites progredes 
along Jupiter’s ecliptic, with pretty nearly the same velocity 
as the node regredcs, assuming these motions to be due to 
the sun’s disturbing force; they are, however, principally 
due to the oblateress of the planet. 

Parallactic Inequality, 

105, In carrying on the approximations to a higher order, 
it is found that the expressions for the forces contain terms 
whose argument is the moon’s elongation^ i,e, the difference 
of longitude of the sun and moon. These terms appear 
in P and in T of the fourth order, and since 1 — m, the 
coefficient of is near unity, the terms will become of the 
third order, and therefore of considerable importance in the 
values of u and of t. We shall work out these expressions 
from the earliest steps ; and we select these terms on account 
of the peculiar use which has been made of them to determine 
the sun’s parallax, whence they have received the name of 
parallactic inequality. 

Let us go back to the expressions for the forces 
Arts, (23), (24). In the first place we may remark that 
the terms we are about to compute n>uat be independent 
of s ; because any term into the composition of which s enters 
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will necessarily have fj in the coefficient of the argument. 
We may therefore at the outset suppose the orbit to coincide 
with the ecliptic. We have then 

SM^ = - ^r'MG cos (0 - &) + MCf\ 

SiP^ V^ / CO8(0 6) ^,3 + ./ ^ cos (0—0 j J 

= 7 ( 1 + + V?' cos2{0-0')}] , 

I 1 r 3/C(r K(7'‘ 1 

«« M» = r'r“’r' {H‘:fcos2(0-^')}J. 

We must now substitute these values in 


D M 

fMG EG\ , , 

/ 1 

1 \ 

P= i/7,-3 +- 

SEV 

Is'J/’ 

SE?) 


— mr 

fcl.- 

-A.'l 

\SM' 

Sb') 


and pick out the terms which have for argument (0 — 0 '). 

y f(MO 7?(/\ , , 

In m 1 J be a term 

3 m co3(0-^)= 7 ^.- cos ( 0 - 0 ). 

In ™'’*’ (^'^3 ~ 003(0 — 0') WO suall have the terms 

' - - CO3(0 - 0') + cos(0- 0') cos2 (0 - 0')}, 


m 


which produce 

and finally mr ~ sin( 0 - 0 ') will produce 

Hi - (i- V)sm(0-0) = i j^^:j:j£Sin(0-0). 



PARALLACTIC INEQUALITY. 


Therefore the new terms will be 

Pu; ^ /ru'^ L + iU ^ ^ A^M ^ 

. T mii^ E-M , K-M , 

in ^...— 2 r ,, , ,Liii (a— t ' )=- 2 m — ,, ,^8111 ( 6 /—^). 

^ IcK' E-\-M ^ ^ ^ n h:^\T V / 


E^M 


Since — (= J ,7 nearly) is of the second order, these terms 
are of the fourth order. 

106. On the principle of the superposition of small 
disturbances, we may compute the ell’ccts of these terms by 
themselves as disturbers of tbe elliptic path. 

[' + - > 

/(V "" ~ J'J+ m7 


2 ^ = a (I cos((^^^- a)+ }, 

— ae sin {cd — a) ; 

therefore 

T clu . . . 


dd ' 


■ 0, to the fout til order, 


- -j- M=a, to the second order, 


f ^ in o ..AJ-Ma 

+ >■ H+Jtz 

Substituting in tbe differential equation for m, we get 

d"u r, ,JS-Ma' ,,, 07. 1 

d7 +“=“ L " Mia 1 • 

Assume « = a |^l + + A ~ cos((l-?>i) 0-/3} +...J > 
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therefore 


therefore u-a 

1+ C0S{(1 m)e 

107. The corresponding term in the value of 6 will also 
bo of the third order. 



1 1 r, 

“a ~ "a 1 • 

u a L 

••••+ J cos{(l- m) 0 - + j , 

C T , 

I /i^u^ is of the fourth order and will not rise in t ; 

therefore 


dt 1 r. 

+ >n)0 , 

<=ir0+ 

F L 

+ +•”] » 

and Q— ft -f . 

- 'im sin {(1 - +••• • 


Mayer was the first who applied this term to the deter- 
mination of the sun’s parallax, by comparing the analytical 
expression of the coefficient with its value as deduced from 

observation. The values of m and of and therefore of 

iP 

^ being pretty accurately known, ~ will be determined, 

that is, the ratio of the sun’s parallax to that of the moon ; 
but the moon’s parallax is well known ; therefore, also, that 
of the sun can be calculated. The value so obtained by 
L-'place was 8*6"; but Hansen, by the discussion of a great 



SECULAR ACCELERATION. 


95 


number of modern observations, has found 8*92", a value 
which agrees very closely with recent determinations by other 
methods. 


Secular Acceleration. 

108. Halley, about 1G93, found, by the comparison of 
ancient and modern eclipses, that the moon’s mean revolution 
is now performed in a shorter time than at the epoch of the 
recorded Chaldean and Babylonian eclipses. This pheno- 
menon, called the secular acceleratmt oj the moon^s mean 
motion^ lias not yet been fully accounted for. For a long 
time it was altogether unexplained, but in 1787 Laplace 
gave what seemed to be a satisfactory and complete account 
of it. 

The value of /?, Art. (54), on which the length of the 
mean period depends, is found, when the approximation is 
carried to a higher order, to contain the quantity e the 
excentricity of the earth’s orbit. Now, this excen tricity is 
undergoing a slow but continual change from the action of 
the planets, and therefore as deduced from observations 
made in different centuries, will have different values. 

The value of i> is at present increasing, or the mean 
motion is being accuiorated, and it will continue thus to 
increase for a period of immense, but not infinite duration ; 
for, as shewn by Lagrange, the actions of the planets on 
the excentricity of the earth’s orbit will be ultimately 
reversed, e' will cease to diminish and begin to increase, and 
consequently p will begin to decrease, and the secular 
acceleration will, so far as this cause is concerned, become a 
secular retardation. 

It is worthy of remark that the action of the planets on 
the moon, thus transmitted through the earth’s orbit, is more 
considerable than their direct action. 
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109. In Ills investigation, Laplace treated the cxccntricity 
e of the earth’s orbit as a constant in the differential equa- 
tions, and only considered its variation in the final result. 
This of course greatly simplified the work, but ought to have 
been looked upon only as a first approximation ; yet, strange 
to say, the result so obtained agreed almost perfectly with 
that which had been deduced from a comparison of ancient 
and modern observations. ‘‘ Malheurcuscment cette iner- 
vcilleusc concordance i\ laquelle Laplace a dfi peut-etre une 
partic de I’dclat do sa plus belle decouverte, est sensiblement 
alt^ree par Ics approximations suivantes.”^* 

Some years ago Professor Adams took up the question 
ab initio without any limitation, and introduced the variability 
of e into the differential equations themselves. This strictly 
correct process gave a value of the secular variation only 
about half of Laplace’s whicli had agreed so well with 
observation. 

'J’hc variation of exccntricity docs not therefore account for 
the whole observed change, and the mean motion of the moon 
is affected by some other cause or causes at present unknown. 
Perhaps the remaining acceleration is only apparent, and 
arises from a gradual retardation of the earth’s diurnal motion. 
The effect would be precisely the sr e, and such a retar- 
dation may be due to the action of the tides. Whatever be 
the cause, the explanation of this unexplained phenomenon is 
a question Avorthy of the mathematician’s most serious efforts. 

Inequalities depend inq on the Figure of the Earth, 

110. The earth, not being a perfect sphere, Avill not at- 
tract as if the whole of its mass were collected at its centre : 
hence, some correction must be introduced to take into 
account this want of sphericity, and some relation must exist 


* roiitecoulaut, dn Momk. Supplement, 18G0, 
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between the oblatcness and the disturbance it produces. 
Laplace in examining its effect found that it satisfactorily 
explained the introduction of a term in the longitude of the 
moon, which Mayer had discovered by observation, and the 
argument of which io the true longitude of the moon’s 
ascending node. 

By a comparison of the observed and theoretical values 
of the coefficient of this term, we may determine the 
oblateiicss of the earth with as great accuracy as by actual 
measures on the surface. 

111. By pursuing his investigations, with reference to 
the oblateness, in the expression for the moon’s latitude, 
Laplace found that it would there give rise to a term in 
which the argument was the true longitude of the moon. 

This term, which was unsuspected before, will also serve 
to determine the earth’s oblatencss, and the agreement with 
the result of the preceding is almost perfect, giving the 
compression which is about a mean between the 

different values obtained by other methods. 

Perturhations dm to Venus, 

112. After the expression for the moon’s longitude had 
been obtained by theory, it was found that there was still 
a slight deviation between her calculated and observed places, 
and Biirg, who discovered it by a discussion of the observa- 
tions of Lahire, Flamsteed, Bradley, and Maskelyne, thought 
it could be represented by an inequality whose period would 
be 184 years and coefficient 15 ". This was entirely con- 
jectural, and though several attempts were made, it was not 
accounted for by theory. 


Pontecoulant, Systme du Mcndtty vol, iv. 


0 
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About 1848, Professor Hansen, of Seeberg, in Gotha, 
having commenced a revision of the Lunar Theory,* found 
two terms, which had hitherto been neglected, due to the 
action of Venus. One of them is direct and arises from 
a ^remarkable numerical relation between the anomalistic 

* motions of the moon and the sidereal motions of Venus 

* and the earth ; the other is an Indirect effect of an Inequality 
‘of long period in the motions of Venus and the earth, 
‘which was discovered some years affo by the Astronomer 
‘lloyal.’^ 

The periods of these two inequalities are extremely long, 
one being 273 and the other 239 years, and their coefficients 
are respectively 27*4" and 23*2". ‘ These arc considerable 

‘ quantities in comparison with some of the inequalities already 
‘recognised in the moon’s motion, and, when applied, they 
‘ are found to account for the chief, indeed the only 
‘remaining, empirical portion of the moon’s motion in 

* longitude of any consequence ; so Aat their discovery may 
‘ be considered as a practical completion of the Lunar Theory, 
‘ at least for the present astronomical age, and as establishing 
‘the entire dominion of the Newtonian Theory and its 
‘ analytical application over that refractory satellite. ’f 

Motion of tlie Ecliptic, 

118. We have seen, Art. (14), that our plane of reference 
is not a fixed plane, but its change of position is so slow 
that we have been able to neglect it, and it is only when 
tne approximation is carried to a higher order that the 
necessity arises for taking account of its motion. 


♦ Report to the Annual (icneral Meeting of the Royal Astronomical Society, 
JFcb. It, 18 IS. 

t Addrc&s of Sir John Hersohel to the Meeting of the Royal Astronomical 
Society. 
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It has been found to have an angular velocity, about an 
axis in its own plane, of 48" in a century, and the correction 
thus introduced produces in the latitude of the moon a term 
— cco cos {6 — <^), 

where ea is the angular velocity of the ecliptic, ^ the angular 

velocity with which the ascending node of the moon’s orbit 
recedes from the instantaneous axis about which the ecliptic 
rotates, (fy the longitude of this axis at time and 0 the 
longitude of the moon at the same instant. 



Let TA771 be the position of the ecliptic at time 
A the point about which It is turning, TA = 

MN the moon’s orbit, 3f the moon, and Mm a perpen- 
dicular to the ecliptic ; 'Tm = 0 ; Mm = lat. = 13. 

{ the inclination of the orbit, and N the longi- 
tude of the node. 

Let APN'in be the ecliptic after a time hi. 

Any point whose longitude is L may be considered as 
moving perpendicularly to the ecliptic with a velocity 
<o sin (L --(f)). 

Hence, the point N will move in the direction PN with 
a velocity a> sin {N—^). And N' will move along PN' with 
a velocity o) sin {N - 4>) cot/; 

dN 

therefore = (o sin {N— (fy) cot/. 

Again, the. point of the ecliptic OO"* in advance of N 
will move towards the moon’s orbit with a velocity 
(o sin {90° + N— (fy) J 
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therefore 


di 

dt 


— « cos {N- <})). 


Now, cot 0 ), and ‘ 


- may be consiuciuu ww- 


stant in integrating ; 

therefore hN ~ cto cos {N— (f>) cot t , 

Si = COD sin (N- 
and if we have 


g'i|r = - 


SN _ 
cos t 


-.cos {N—<}>), 


Now, sin ^ == sin i. sin yfr • 
cos )8 . S/3 = cos t . sin + sin i. cos . S>|r 

= 00 } {cosisIn-^/r sin(i\r~ — ^o^yjr cqs{N — 4 >)]} 
but cosi.sin'^ =cosy8 8in(0— and cos^ = cos^ cos(0—N^) ; 
therefore S/3 = — c(d cos (0 — 

The discovery of this term is due to Professor Hansen; 
its coefBcient is extremely small, about 1*5"; but, being of 
a totally different nature from those due to successive 
approximations, it was thought desirable to examine it, and 
the above investigation, which was communicated to me by 
J. C. Adams, Esq.,* will be read with interest on account 
of its elegance. 

With respect to tho foregoing investigation, perhaps the following remarks 

(ly 

will not bo superfluous; — Tho value of is not the actual velocity of W, 

but its velocity relatively to the position of the node as determined when tho 
motion of the ecliptic is neglected ; its integral is therefoi*e cN the change of 
longitude due to this motion, and in this integration no constant is added, zero 

(LV 

being taken for the mean value. The periodic forms of both and oN shew 
that they oscillate about mean values, the time of a complete oscillation being 


* Now Lowndcan Piofcssor of Astronomy in the University of Cambridge. 
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that required by sin {N — tp) and cos (X—<p) to go througli their cycle. This 
relative motion of the node is analogous to that of a particle moving in a straight 
line under the action of a force varying directly iis the distance. 

Similar remarks apply to the inclination. 

It must filso be borne in mind that the result obtained ^/3 is not the difference 
between the latitude referred to the actual ecliptic and that referred to a fixed 
plane, but is the difference Ixitwcen the calculated latituile referred to the actual 
ccliiDtic on supposition of its being fixed, and the correct latitude referred to the 
same actual ecliptic when its ra'-tion is taken into account. 


We may obtain an approximate value of tbe coefficient 

cco by substituting for it , where a is the number of seconds 

through which the ecliptic is deflected in one year = 0*48", 
and n is the number of years in which the node of the moon’s 

orbit makes a complete revolution = 18 0 ; for then, is 
the angle described by the node in one year; therefore 
IS the ratio of ft) : ~ ~ , supposing ^ to remain con- 

stant, which is nearly the case ; therefore 


Cft) = 


na 

2ir 


9*3 X 0*48" 
3*14 


= 1*42".«* 


* This affords the solution of a problem proposed in the SenatC'Houso 
in the January Examination v. Question 21, Jan, 22. 
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CHAPTER VIIL 

HISTORY OF THE LUNAR PROBLEM BEFORE NEWTON. 

114. The idea wliich most probably suggested itself to 
the minds of those men who first considered the motion of 
the moon among the stars, was that it described a circle with 
uniform velocity about the earth as a centre ; and this first 
rough result is represented in our equations by neglecting 
all small terms and writing Q ==p^, w = a. 

It must, however, have been very soon perceived that 
the actual motion is far from being so simple, and that 
the moon moves with very difierent velocities at different 
times. 

115. The earliest recorded attempts to take into account 
the irregularitcs of the moon’s motion were made by 
Hipparchus (140 b.c.). lie imagined the moon to move with 
uniform velocity in a circle, of which the earth occupied, 
not the centre, but a point nearer to one side. By a similar 
hypothesis ho had accounted for the irregularities in the 
sun’s motion, and his success in this led him to apply it also 
to the moon. 

It is clear that, on this supposition, the moon would seem 
to move faster when nearest the earth or in perigee, and 
slower when in apogee, than at any other points of her 
orbit, and thus an apparent unequal motion would be 
produced. 
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Let BAM be a circle, CA a radius, E a point in A 0 
near (7; (7 j5, ED two parallel 
lines making an angle a with 
CA. 

Suppose a body M to describe 
this circle uniformly with an 
angular velocity jo, the time being 
reckoned from the instant when 
the body was at 5, and the longi- 
tude as seen from E being reckoned 
from the line ED ; 
therefore DEM = BCM - 

AEM = e-oL, ACM^pt--a. 

EC . 

Now is a small fraction, and if wo represent it by 
we shall have 

sinilf=55;8in AEM 
CM 

= e sin (0 — a) 

= e sin + M- a), 

^ csiu(;;^-a) 

tanJf= - -E-; 

1 —ecos (jft’-a) 

this would give M, an^^ then 0 by the formula 0—pt-\rM. 

This was called an excentric^ and the value of e was called 
the excentricitij^ which, for the moon, llipparchus fixed at 
sin 5 1 . 


116. Another method of considering the motion was by 
means of an epicycle.^ which led to the same result. 

A small circle PJ/, with a radius equal to EC of previous 
figure, has its centre in the circumference of the circle RED 
(which has the same radius as that of the exccntric), and 
moves round E with the uniform angular velocity p^ the 
body M being carried in the circumference of the smaller 
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circle, the radius PM remaining parallel to itself, or, which 
is the same thing, revolving 
from the radius PE with the 
same angular velocity jp, so 
that the angle EPM equals 
PEA. 

Now, when the angle AEP 
equals the angle ACM of the 
former figure, it is easily seen 
that the two triangles EPM^ 

ECM arc equal, and there- 
fore tlic distance EM and the angle AEM will be the same 
in both, and the two motions will be identical. 

117. The value of e being small, we find, rejecting &c., 
il/==e sin — a), 

therefore e=:^pt^e sin {pt - a). 

If we reject terms of the second order in our expression 
for the longitude, and make c = 1, we get, Art. (55), 

0 =zpt -f- 2e sin {pt - a), 

which will be identical with the above if we suppose the 
cxccntricity of the excentric to be double that of the elliptic 
orbit. 

Ptolemy (a.d. 140) calculated the excentricity of the 
moon’s orbit, and found for it the same value as Hipparchus, 
viz. 

sin 5° 1' = -i^ 2 j i^carly. 

The excentricity in the elliptic orbit is, we know, about 
and ^ and will pretty nearly reconcile the two values of 6 
given above. This shews us, that for a few revolutions the 
moon may be considered as moving in an exccntric, and her 
positions in longitude calculated on this supposition will be 
correct to the first order. 
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Her distances from the earth will not however agree ; for 
the ratio of the calculated greatest and least distances would 

be or } f, while that of the true ones would be ^ 


It would, therefore, have required two different cxcentrics 
to account for the changes in the moon’s longitude and in 
her radius vector. Changes in the latter could not, however, 
be easily observed with the rude instruments the ancients 
possessed, and it was very long before this inconsistency was 
detected. 


118. We have said that the moon’s longitude, calculated 
on the hypothesis of an exceutric, will bo pretty accurate for 
a few revolutions. 

The data requisite for this calculation are, the mean 
angular motion of the moon, the position of the apogee, and 
the magnitude of tlie exccntricity. 

But it was known to Hipparchus and to the astronomers of 
bis time, that the point of the moon’s orbit where she seems 
to move slowest, is constantly changing its position among 
the stars. Now this point is the apogee of llipparchus’a ex- 
centric, and he found that he could very conveniently take 
account of this further change by supposing the excentric 
itself to have an angular motion about the earth in the same 
direction as the moon, so as to make a complete revolution 
in about nine years, or about 3° in each revolution.^ 

This motion of the apsidal line follows also from our 
expression for the longitude, as shewn in Art. (75). It is 
there, however, connected with an ellipse instead of an 


* On tlio supposition of an cjncyrle, this motion of the aj>so coulfl as e.'isily 
be represented by supposing the radius which connects tfie moon with the rjentro 
of tlie ef)Icycic to have this unih»rm angular velocity ‘-f about ‘o' in each revolu- 
tion, and cil-o in the '^alne <1imcliou. i 
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cxccntric ; and though the discovery that the clUptlc ns the 
true form of the fundamental orbit was not the next in the 
order of time after those of Hipparchus, yet, as all the irre- 
gularities which were discovered in the intervening seventeen 
centuries are common both to IIippa*'chus’s cxccntric and to 
Kepler’s ellipse, It will be as well for us to consider at once 
this new form of the orbit. 

Elliptic Form of the OrhiL 

119. We need not dwell on the steps which led to this 
great and important discovery. Kepler, finding that tlie 
predicted places of tlie planet Mars, as given by the circular 
theories then in use, did not always agree with the computed 
ones, sought to reconcile these variances by other combina- 
tioiis of circular orbits, and after a great number of attempts 
and failures, and eight years of patient investigation, he 
found it necessary to discard the cxccntrics and epicycles 
altogether, and to adopt some new supposition. An ellipse 
with the sun in the focus was at last his fortunate hypothesis? 
which was found to give results in accordance with obser- 
vation; and this form of the orbit was, with equal success, 
afterwards extended to the moon: buw the departures from 
elliptic motion, due to the disturbing force of the sun, arc, 
in the case of the moon, much greater than the disturbances 
of the planet Mars by the other planets. 

In Kcplers’s hypothesis, then, the earth is to be considered 
as occupying the focus of an ellipse, in the perimeter of which 
the moon is moving, no longer with either uniform linear or 
angular velocity, but in such a manner that the radius vector 
sweeps over equal areas in equal times. 

This agrees with our investigation of the motion of two 
bodies, Art. (10). 
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Evection, 

120. The hypothesis of an excentric, whose apse line has 
a progressive motion, as conceived by Hipparchus, served to 
calculate with conside^'able accuracy the circumstances ot 
eclipses; and observations of eclipses, requiring no instru- 
ments, were then the only ones which could be made with 
sufficient exactness to test the truth or fallacy of the sup- 
position. 

Ptolemy (a.d. 140) having constructed an instrument, by 
means of which the positions of the moon could be observed 
in other parts of her orbit, found that they sometimes agreed, 
but were more frequently at variance with the calculated 
places; the greatest amount of error always taking place at 
quadrature and vanishing altogether at syzygy. 

What must, however, have been a source of great per- 
plexity to Ptolemy, when he attempted to investigate the 
law of this new irregularity, was to tind that it did not 
return in every quadrature, — in some quadratures it totally 
disappeared, and in others amounted to 2^* 39', which was 
its maximum value. 

By dint of careful comparison of observations, he found 
that the value of tins second inequality in quadrature was 
always proportional to that of the first in the same place, 
and was additive or subtractive according as the first was 
so ; and thus, when the first inequality in quadrature was at 
its maximum or 5° 1', the second increased it to V 40', which 
was the case when the apse line happened to be in syzygy 
at the same time.* 


* It would seem as if Hipparchus had felt the necessity for some further 
modification of lus first hypothesis, though he was unable to determine it; for 
there is an observation made by him on the moon in the position here specified 
when the error of his tables would be greatest ; and at a time also when she 
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But if the apse line was in quadrature at the same time as 
the moonj the second iiuiquality vanished as well as the first. 

The mean value of the two inequalities combined was 
therefore fixed at G*" 20;^, 

121. To represent tliis new inequality, which was sub- 
sequently called the Ptolemy Imagined an cxcentric 

in the circumference of which the centre of an epicycle moved 
while the moon moved in the circumference of the epicycle. 

The centre of the cxcentric and of the epicycle he supposed 
in syzygy at the same time, and both on the same side of 
the earth. 



Thus, if E represent the earth, 


S sun, 

M moon, 


0 the centre of the cxcentric UKT in syzygy, 
jB, the centre of the epicycle, would also be in 
syzygy. 

Now conceive c, the centre of the cxcentric, to describe a 
small circle about E in a retrograde direction cc', while iZ, 
the centre of the epicycle, moves in the opposite direction, 


was in tlic nonngesimaJ, so that any error of longitude, arising from her yet 
uncertain parallax, would be avoided. Ptolemy, who records the observation, 
employs it to calculate the evoction, and obtains a result agreeing with that 
of his own obsen’atious. (Soc Dclambre. .1.^. 



Ptolemy’s hypothesis. 


lOJ) 

in such a manner that each of the angles S'Ec^ S'ER may 
be equal to the synodical motion of the moon, that is, her 
mean angular motion from the sun ; SES' being the motion 
of the sun in the same time. 

Now we have seen, Art. (llC), that the first inequality 
was accounted for by supposing the epicycle liM to move 
into the position rm^ r and 11 being at the same distance 
from E^ and rm parallel to iiJ/,* the first inecpiality being 
the angle rEm. But when the centre of the epicycle is at 
and RM' is parallel to rm^ the inequality becomes R EM' ^ 
and we have a second cori ection or inequality inEM', 

122. That this hypothesis will account for the phenomena 
observed by Ptolemy, will be readily understood. 

At syzygics, whether conjunction or opposition, the centres 
of the exceiitric and epicycle arc in one line with the earth 
and on the same side of it ; the points r and R coincide, as 
also m and M'. Hence inEM' 

At quadratures (figs. 1 and 2) d and R arc in a straight 
line on opposite sides of the earth, and therefore R and r at 

riff. 2. 




their furthest distance. If, however, M' and m be at the 
Bame time in this line, or, in other words, if the apse line 


* For simplicity we leave out of consideratir n the motion of the apse. 
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be in quadratures (fig. 1), the angle mEM' will still be zero, 
or there will be no error in the longitude. But, if the apse 
line is in syzygy (fig. 2), the angle mEM' attains its greatest 
value.* 

Ptolemy, as we have said, found this greatest value to be 
2^ 39', the angle mEr being then 5° 1'. 

123. Copernicus (a.d. lo 13), having seen that Ptolemy’s 
hypothesis gave distances totally at variance with the obser- 
vations on the changes of apparent diameter,! made another 
and a simpler one which accounted equally well for the in- 
equality in longtitudc,"and was at the same time more correct 
in its representation of the distances. 

Let E be the earth, OD an epicycle whose centre C de- 
scribes the circle O' CO" about E with the moon’s mean 
angular velocity. 



Let CO, a radius of this epicycle, be parallel to the apse 


♦ If Ptolemy had used the hypothesis of an excentric instead of an ejncycle 
for the first inc'quulity of the moon, an epic\cle would have reprtsented the 
second inequality more simply than his method did. Dr. Whewells History of the 
Inductive ^Sciences, vol. I., p. 230. 

f See Delambre, Moderncy vol. I., p. 116. WhewelVa History of Inductive 
Sciences^ vol. i., p. 305. 
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Ill 


line EA^ and about 0 as centre let a second small epicycle 
be described, the radii CO and 021 ho so taken that 


CO - OM 
CE 


= sin 5“ 


r, and 


CO+ OM 
CE 


= sin7° 40’. 


The radius OM muse now be made to revoke from the 
radius OC twice as rapidly as EC moves from E8^ so that 
the angle COM may be always double of the angle CE8. 

From this construction, it follows that in syzygies the 
angle CE8 being 0° or 180°, the angle C021 is if or 3G0°; 
and therefore C and M arc at their nearest distances, as in 
the positions C" and (7"' in the figure. Then CM= CO - OM^ 
and the angle CE2I will range between 0° and 5° 1', the 
greatest value being attained when the apse line is in quad- 
rature. 

When the moon is in quadrature CES—^if or 270°, 
and, therefore, COM=\^if or 540° and C and M are at 
their greatest distance apart, as in the position (7"; then, 
(7iT/= 6^ 0 + Od/j and ihe angle CEM will range between 
0° and T 40', the former value when the apse line is itself 
in quadrature, and the latter when it is in syzygy. 


124. Thus the results attained by Ptolemy’s construction 
are, as far as the longitudes at syzygies and quadratures are 
concerned, as well represented by that of Copernicus; and 
the variations in the distances of the moon will be far more 
exact, the least apparent diameter being 28' 45" and the 
greatest 37' 33" ; whereas, Ptolemy’s would make the greatest 
diameter T.* 

The values which modern observations give vary between 
28' 48" and 33' 32". 


♦ Delambre, M<>dtrne. 
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125. It will not now be difficult to shew that the intro'- 
duction of tliis small epicycle corresponds with that of the 
term sin {(2 — 2ni - c) - 2/Q -f a} in our value of 0. 

For, referring to the preceding figure, we have 

OEM^ sin OEU - sin OME 

OE 

\J Tj 

= -/ sin (2 . SEC - A EM) 

- (moon’s mean long. - sun’s long.) 

— (moon’s true long. ~ long, of apse)}, 
and OEM being a small angle whose maximum is T 
we may write moon’s moan longitude ihstoad of the true 
in the argument, and also KC for therefore, 

OM 

sin (2 (moon’s mean longitude - sun’s longitude) 

— (moon’s mean longitude — longitude of apse)} 
= 70^' sin [2 [pi - ( 7 /?^;^ + /3)j - (I - a}j 
= 4770" sin [(2 — 2;/^ — c)pt — 2/3 4- a}. 

The value of the coefficient is from modern observations 
found to be 4589*61". 

120. In Art. (77), we have considered the effect of this 
second inequality in another light, not simply as a small 
quantity additional to the first or elliptic inequality, but as 
forming a part of this first; and therefore, modifying and 
constantly altering the cxcentricity and the uniform pro- 
gression of the apse line. 

BouHiaud (A. D. 1645), by whom the term Evection was 
first applied to the ccond inequality, seems to hint at some- 
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lia 

thing of this kind in the rather obscure explanations of Jus 
lunar hypothesis, which, never having been accepted, it would 
be useless to give an account of it.* 

In Ptolemy’s theory. Art. (121), the evectlon was the 
result of an apparent increase of the first lunar epicycle 
caused by its approaching the earth at quadratures ; but, in 
this second method, it is the result of an actual change in 
the elements of the elliptic orbit. 

D’Arzachel, an Arabian astronomer, who observed in 
Spain about the year 1080, seems to have discovered the 
unequal motion of the apsides, but his discovery must have 
been lost sight of, for Ilorrocks, about 1G40, re-disco vered it 
‘ in consequence of his attentive observations of the lunar 

* diameter : he found that when the distance of the sun from 

* the moon’s apogee was about 45° or 225°, the apogee was 
‘ more advanced by 25° than when that distance was about 
‘ 135° or 315°. The apsides, therefore, of the moon’s orbit 
‘ were sometimes progressive and sometimes regressive, and 
^required an equation of 12° 3(f, sometimes additive. to their 
^ mean place and sometimes subtractive from it.’f 

Ilorrocks also made the cxccntrlcity variable between the 
limits *06686 and *04362. ^ 

The combination of these two suppositions was a means 
of avoiding the introduction of Ptolemy’s exccntric or the 
second epicycle of Copernicus : their joint effect constitutes 
the evection. 


♦ AprSs avoir etabli lea mouvements et les epoquoa de la lane, Boulliaud revient 
^ I’explication de Icvection ou de la scconde incgalite. Si ea theorie u’a pas fait 
fortune, le nom du moins cst reste. ‘ Kn meine temps quo la lune avance sur son 
‘ cone autour de la terre, tout le sj/steme de la lune ent dt'place j la terre omj)ortant la 
‘lune, rejette loin d’elle Tapogee, et rapproche d'autant le perigee; mais cette 
‘ Evection d. des homes fixees.’ Delambre, Jlist. <h V Ast, Mod., tom. ii. p. 157. 
t Small’s Astronomical Discoveries of Kepler, p. 307. 


Q 
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Variation. 

127. After the discovery of the evection by Ptoleray, a 
period of fourteen centuries elapsed before any further ad- 
dition was made to our knowledge of the moon’s motions. 
Hipparchus’s hypothesis was found sufficient for eclipses, and 
when corrected by Ptolemy’s discovery, the agreement be- 
tween the calculated and observed places was found to 
extend also to quadratures; any slight discrepancy being 
attributed to errors of observation or to the imperfection 
of instruments. 

But when Tycho Brah4 (a. d. 1580) with superior instru- 
ments extended the range of his observations to all inter- 
mediate points, he found that another inequality manifested 
itself. Having computed the places of the moon for differ- 
ent parts of her orbit and compared them with observation, 
he perceived that she was always in advance of her com- 
puted place from syzygy to quadrature, and behind it from 
quadrature to syzygy ; the maximum of this variation taking 
place in the octants, that is, in the points equally distant 
from syzygy and quadrature. The moon’s velocity therefore, 
so ^ar as this inequality was concerned, was greatest at new 
and full moon, and least at the first and third quarters.* 


* ‘It appears that Mohammed- Aboiil-Wefa-al-Boiirdjani, an Arabian astro- 
nomer of the tenth century, who resided at Cairo, and observed at Bagdad in 975, 
‘ discovered a third inequality of the moon, in addition to the two expounded by 
‘ Ptoleray, the equation of the centre and the evection. This third inequality, the 
‘variation, is usually supposed to have been discovered by Tycho Brahe, six 

‘centures later In an almagest of Aboul-Wefa, a part of which exists in 

‘ the Royal Ijibrar}’^ at Paris, after describing the two inequalities of the moon, 
‘he has a Section IX., “Of the third anomaly of the moon called Muhnznt or 

^ Prosnensis"" But tliis discovery of Aboul-Wefa appears to have excited 

‘ no notice among his contemporaries and followers ; at least it had been long 
‘ quite forgotten, when Tycho Brahe re-discovered the same lunar inequality.’ 
Whewell's Hist, of Induct ive /^••ienreSf vol. i. p. 243, 
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Tycho fixed the maximum of this inequality at 40' 30". 
The value which results from modern observations is 
39 ' 30 ". 


128 . We have already two epicycles, or one epicycle and 
an excentric, to explain the first two inequalities: by the 
introduction of another epicycle or excentric, the variation 
also might have been brought into the system ; but Tycho 
adopted a ditFcrent method:* like Ptolemy, he employed 
an excentric for the evection, but for the first or elliptic 
inequality be employed a couple of epicycles, and this 
complicated combination, which it is needless further to 
describe, represented the change of distance better than 
Ptolemy’s. 

To introduce the variation^ he imagined the centre of the 
larger epicycle to librate backwards and forwards on the ex- 
centric, to an extent of 40^^' on each side of its mean position ; 
this mean place itself advancing uniformly along the excentric 
with the moon’s mean motion in anomaly ; and the libration 
was so adjusted, that the moon was in her mean place at 
syzygy and quadrature, and at her furthest distance from it 
in the octants, the period of a complete libration being half 
a synodical revolution. 

Annual Equation. 

129 . Tycho Brahd was also the discoverer of the fourth 
inequality, called the annual equation. This was connected 
with the anomalistic motion of the sun, and did not, like the 
previous inequalities, depend on the position of the moon ’U 
her orbit. 


♦ For a full description of Tycho’s hypothesis, see Delambre, Hiixt. de V Ast, 
Mod. tom. I. p. 102, and An Account the, Astronomical Discoveries of Kepler 
by Robert Small, p. 130. 
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Having calculated the position of the moon corresponding 
to any given time, he found that the observed place was 
behind her computed one while the sun moved from perigee 
to apogee, and before it in the other half year. 

Tycho did not state this distinctly, but he made a cor- 
rection which, though wrong in quantity and applied in an 
indirect manner, shewed that he had seen the necessity and 
understood the law of this inequality. 

He did not try to represent it by any new cxcentric or 
epicycle, but he increased by (8m. 13s.) 9\n[siin^s anomaly) 
the time which had served to calculate the moon’s place 
thus assuming that the true place, after that interval, would 
agree with the calculated one. Now, as the moon moves 
through 4' 30" in 8m. I3s., it is clear that adding (8m. ISs.) 
m\[sun''s a^iomnly) to the time is the same thing as sub- 
tracting (4' 30") sin (.9«w’5 anomaly) from the calculated lon- 
gitude, which was therefore the correction virtually intro- 
duced by Tycho.f Modern observations shew the coefEcIent 
to be 11' 0". 

We have seen. Art. (82), how this inequality may be 
inferred from our equations. 

Beduction, 

130. The next inequality in longitude which we have to 
consider is not an inequality in the same sense as the fore- 
going r that Is, it does not arise from any irregularity in the 
motion of the moon herself in her orbit, but simply because 
that orbit is not in the same plane as that in which the 


* That is, tlie equation of time which he used for the moon differed by thai 
quantity from that used for the sun. 

t Uorrocks made the correction in the same manner as Tycho, but sc 

increased it that the coiTcsponding coefficient was lU 51" instead of 4' ilO" 
Flamsteed was the first to r pply the correction to the longtitude instead of the time 



REDUCTION, 


117 


longitudes are reckoned, so that even a regular motion in 
the one would be necessarily irregular when referred to the 
other. Thus if NMn be the moon’s orbit and yNm the 
ecliptic, and if M the moon be re- 
ferred to the ecliptic by the great 
circle Min perpendicular to it, then 
3IN and viN are OO*", 180°, 270°, 
and 3G0° simultaneously, but they 
differ for all intermediate values: the diff’ercnco between 
them is called the reduction. 

The difference between the longitude of the node and that 
of the moon in her orbit being known, that is the side NM of 
the right-angled s{)herical triangle NMiiij and also the angle 
N the inclination of the' two orbits, the side Nrn may be 
calculated by tlnj rules of spherical trigonometry, and the 
difference between it and NM^ applied with a proper sign 
to the longitude in the orbit, gives the longitude in the 
ecliptic. ^ , 

Tycho was the first to make a table of the reduction 
instead of calculating the spherical triangle, llis formula 
was 

reduction = sin 2 Zf — ^ tan*^/ sln4Zf, 

where / is the Incline ion of the orbit and L the longitude of 
the moon diminished by that of the node. 

The first term corresponds with the term — sin2 {gpt — 7 ) 
of the expression for 6, 

Latitude of the Moon. 

131. That the moon’s orbit is inclined to the ecliptic was 
known to the earliest astronomers, from the non-recurrence 
of eclipses at every new and full moon ; and it was also 
known, sinces the eclipses did not always take place in the 
same parts of the heavens, that the line lof nodes represented 
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by Nn^ in the preceding figure, has a retrograde motion on 
the ecliptic, N moving towards T. 

Hipparchus fixed the inclination of the moon’s orbit to 
the ecliptic at 5®, which value he obtained by observing the 
greatest distance at which she passes to the north or south 
of some star known to be in or very near the ecliptic, as for 
instance the bright star llegulus ; and by comparing the re- 
corded eclipses from the times of the Chaldean astronomers 
down to his own, he found that the line of nodes goes round 
the ecliptic in a retrograde direction in about 18| years. 

Tliis result is indicated in our expression for the value 
of the latitude by the term k sin (gd — 7), as we have shewn 
Art. (85). 

182. Tycho Brahd further discovered that the inclination 
of tlie lunar orbit to the ecliptic was not a constant quantity 
of 5° as Hipparchus had supposed, but that it bad a mean 
value of 5° 8', and ranged through 9' 30" on each side of 
this, the least inclination 4° 58^' occurring when the node 
was in quadrature, and the greatest 5° 17^' being attained 
when the node was in syzygy.* 

He also found that the retrograde motion of the node was 
not uniform; the mean and the true T^osltions agreed very 


* Ebn Jounis, an Arabian astronomer (died a.d. 1008), whose works were 
translated about 50 years since by Mens. Sedillot, states that the inclination of the 
moon’s orl'it had been often observed by Aboul-Hassan-Aly-ben-Amajour about 
the year 918, and that the results he had obtained were generally greater than 
the 6° of Hipparchus, but that they varied considerably, 

Ebn Jounis adds, however, that he himself had observed the inclination 
seveial times and found it 5° 3', which leads us to infer that he always observed 
ill similar cii'curastanceg, for otherwise a variation of nearly *23' could scarcely 
have escaped him. Sec Delaiubre, Hist, de V Ast. du Moyen Aye, p. 139. 

The mean value of the inclination is 5^ 8' 55*40", — the extreme values are 
4° 67' *22" and 5° 20' 6". 

The me.-in daily motion of the line of nodes is 3' 10*64", or one revolution in 
6793*29 days, or 18 y. ‘218 d. 21 h. 2‘2m. 46 s. 
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well when they were in syzygy or quadrature, but they 
were T 46' apart in the octants. 

By referring to Art. (87) we shall sec that these correc- 
tions, introduced by Tycho Brahe, correspond to the second 
term of our expression for s. 

Since Hipparchus could observe the moon with accuracy 
only in the eclipses, at which time the node is in or near' 
syzygy, we see why he was unable to detect the want of 
uniformity in the motion of the node. 

133. To represent these changes in the position of the 
moon’s orbit, Tycho made the following hypothesis. 

Let JiNF be the ecliptic, K its pole, BA C a small circle, 



having also K for pole and at a distance from it equal 
to S’’ 8'. Then, if we suppose A the pole of the moon’s 
orbit to move uniformly in the small circle and in the 
direction BAG^ the node which is at 90® from both A 
and AT, will retrograde uniformly on the ecliptic, and the 
inclination of the two orbits will be constant and equal to AK, 
But instead of supposing the pole of the moon's orbit to 
be at A^ let a small circle cibed be described with A as polo 
and a radius of 9' 30" ; and suppose the pole of the moon’s 
orbit to describe this small circle with double the velocity 
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of the node In Its synodical revolution which is acconi{)lishcd 
in about 31G days, in such a manner that when the node is 
in quadrature the pole may be at a, the nearest point to 
and at c the most distant point when the node comes to 
syzygy, at d in the first and third c jtants, and at h in the 
second and fourth, so as to describe the small circle in about 
173 days, tlic centre A of the small circle retrograding mean- 
while with its uniform motion. 

By this method of representing the motion, we sec that 


when node is in quadrature' 

«jzygy 


the inclina- 
tion of the 
orbit is 




while at the octants it has its mean value 


= = = 8 '. 

Again, with respect to the motion of the node, since N is 
the pole of KaAc^ it follows that when in syzygy and quad- 
rature, the node occupies its mean place ; in the first and 
third octants, the pole being at the node is before its 
mean place by the angle c/AT^l = (9' 30") coscco’ 8' = 1° 46', 
nearly, and it is as much behind in the second and fourth 
octaLcS. 

So that the whole motion of the node, and the correction 
which Tycho had discovered, were properly represented by this 
hypothesis, which is exactly similar to that which Copernicus 
had im..giucd to explain the precession of the Equinox. 
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SELECTION OF EXAMINATION QUESTIONS 

FROM COLLEGE AND SENATE-IfOUSE PAPERS AND FROM THE 
MODERATORSHIP \ND FELLOWSHIP EXAMINATIONS AT 
TRINITY COLLEGE, DUBLIN. 

1. Define the plane of the ecliptic and prove that, as 

seen from the earth, sin = sinV nearly; \ being 

the latitude of the sun and that of the moon. 

2. Obtain the differential equation of the moon’s radius 
vector. 


3. What is the principle on which the successive approxi- 
mations to the nw)on’s path are obtained ? 

Compare the coeflScients of the principal term in the cen- 
tral disturbing force with the principal term in the central force. 


4. 


If in solving thd equation 


fFu _ P 
dd^'^^~hV 


T du 
h'u’ dd 



we take as our first approximation u — a {l+c co8(^-a)}, 
we obtain at our higher approximations terms containing 0 in 
their coefficient, Kow is this defect avoided by taking 
u = a {1 cos(c^ — a)} as our first approximation? 

Upon what principle do we approximate to the value of o? 


T 

5. Shew that in the expansion of there will be a 

term — sin {(2 — 2m — 2c) 9 — 20-^ 2a} which w«uld 

rise to the second order in the longitude. 

Shew further that the terra with the same argument does 
not appear actually in the longitude — assuming 

w =a [1 + e cos(c^ - a) +...+ ^me cos {(2 - 2m — c) ^ — 2/3 -f a} 
+...4- cos {(2 ~ 2m- 2c) 0 — 2/34-2a}4*...]. 
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G. The longitude of the moon contains a term of the form 
A sin {(2 — 2m — c) yt “2/34- a} 

where A is of the second order; find how it must have 
arisen a..d determine A. 

7. Investigate the following expression for the moon’s 
longitude as far as the second order, supposing the orbits 
of the moon relative to the earth, and of the earth relative 
to the suiij to be originally circles in the plane of the ecliptic, 

6 'pi y 'm? sill {(2 — 2m) pi — 2/3}. 

8. Suppose that whatever constants are involved in the 
arguments 0, </>, yjr^ &c. in the equation 

Lr= A + B sin 0 + G sin <!> + J) sin i/r 4- &c. 
arc perfectly known, and that observation gives the numerical 
values of a considerable number of longitudes Zj, vSre., 

corresponding to the known angles &c., &c., 

v^’^c., &c., shew how the constants A, 7>, C\ &c. 
may be numerically determined in the two distinct cases 
where the observations at our disposal are and are not 
unlimited in number. 

0. Ex; lain how this process fails and may be modified 
when two of the terms in question B sin 6y O sin (f> arc nearly 
synchronous in their periods. 

10. Shew that the evection in longitude, viz. 

y sin {(2 — 2/n — c) — 2/3 + aj, 
may be represented as the joint effect of certain periodic 
changes in the exccntricity of the lunar orbit and in the 
mean longitude of its apse. 

11. Assuming the usual notation of the Lunar Theory, 
explain the physical meaning of the following equations : 

s — k sin {0 — 7), 
vS* = k sin ((jO — 7), 

s = k sin (qd — 7) 4 ;s nik sin {(2 — 2m — g) ^ — 2/347}. 

12. Explain the effect of the term ir^a cos l(2~2;n) ^-“2/31} 
in the moon’s radim vector; find the number of days in 
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the period of the rosiilting inequality and the ratio of the 
axes of the oval orbit. 

13. Investigate the effect of the annual equation on the 
length of the lunar month. 

Assuming r— *017, sin l.T30'=*233 lo, and the mean 

period of a sidereal revolution of the moon = 27d. 7h. llUn, 12s., 
find the difference between a winter and a summer month. 

14. Shew that the moon’s orbit in space is always concave 
towards the sun. 

Assuming that tlie earth describes a circle about tlie sun, 
and the moon a cinde about the cartli and in the same plane 
with the earth’s orbit, compare the curvatures of the moon’s 
orbit in space at her perilielion and aphelion. 

15. Prove that if we go to the third order of approxi- 

mation, the motion of the moon’s apse in one revolution 
of the moon equalj (1 — 3(>0.° 

P T 

10. The expressions for and ,, contain terms of 

the form 3^1 cos and A sin(0— 0') respectively; com- 
pute the terms resulting from them in the values of ti and 0, 

! 

17. Copernicus represented the evectlon by an epicycle 
supcrpos(?d on the epicycle which represented the elliptic 
inequality. Prom the construction given by him fihew 
geometrically that the evection vanishes when the moon’s 
mean elongation from the sun is half her true anomaly, 
and, in general, varies nearly as the sine of twice the 
difference of these angles. 

18. The secular change of inclination of the actual ecliptic 
to a fixed ecliptic being a" annually, show that this will 
give rise to an inequality in the moon’s latitude whose typo 

Is — a" cos 0, nearly; ii being the number of years in 

which the moon’s nodes make a complete revolution, and 
the difference of longitude of the moon and of the ascending 
node of the actual ecliptic. 

W. MKK \LJ-i: AMi 1M:IATI:US CVMIUUJX.li 
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